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Abstract 

We conjecture that all CP violations (both Dirac and Majorana types) arise from a common origin 
in neutrino seesaw. With this conceptually attractive and simple conjecture, we deduce that \x—t 
breaking shares the common origin with all CP violations. We study the common origin of fi—r and 
CP breaking in the Dirac mass matrix of seesaw Lagrangian (with right-handed neutrinos being 
[i—T blind), which uniquely leads to inverted mass-ordering of light neutrinos. We then predict a 
very different correlation between the two small fi—r breaking observables 13 — 0° and # 23 — 45° , 
which can saturate the present experimental upper limit on 8 13 . This will be tested against our 
previous normal mass-ordering scheme by the on-going oscillation experiments. We also analyze 
the correlations of 9 13 with Jarlskog invariant and neutrinoless /3/3-decay observable. From the 
common origin of CP and fi~r breaking in the neutrino seesaw, we establish a direct link between 
the low energy CP violations and the cosmological CP violation for baryon asymmetry. With 
these we further predict a lower bound on #13 , supporting the on-going probes of #13 at Daya 
Bay, Double Chooz and RENO experiments. Finally, we analyze the general model-independent 
Z2 €5 Z 2 symmetry structure of the light neutrino sector, and map it into the seesaw sector, 
where one of the corresponds to the fj,— r symmetry Z^ 1 " and another the hidden symmetry 
Z| (revealed in our previous work) which dictates the solar mixing angle 9\i . We derive the 
physical consequences of this Z| and its possible partial violation in the presence of fi— r breaking 
(without or with neutrino seesaw), regarding the 6\2 determination and the correlation between 
/.£ — t breaking observables. 
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1. Introduction 



We conjecture that all CP violations (both Dirac and Majorana types) arise from a common origin 
in neutrino seesaw. With this conceptually attractive and simple conjecture, we deduce that [i — r 
breaking shares the common origin with all CP violations, since the ll — t symmetric limit enforces 
vanishing mixing angle #13 and thus Dirac CP conservation. 

In a recent work [I], we studied the common origin of soft /U— r and CP breaking in the neutrino 
seesaw, which is uniquely formulated in the dimension-3 Majorana mass term of singlet right-handed 
neutrinos. This formulation predicts the normal mass ordering (NMO) for light neutrinos. In this 
work, we study in parallel a different realization of the common origin of fx — r and CP breaking in 
the "jU— t blind seesaw", where the right-handed neutrinos are singlet under the //— r transformation. 
We then find the Dirac mass-matrix to be the unique place for the common origin of [i — t and CP 
breaking in the fj,—T blind seesaw. Since the Dirac mass-matrix arises from Yukawa interactions with 
Higgs boson(s), this can also provide an interesting possibility of realizing spontaneous CP violation 
with CP phases originating from the vacuum expectation values of Higgs fields. Different from our 
previous construction [T] , we reveal that the common origin of \i — r and CP breaking in the Dirac 
mass- matrix uniquely leads to the inverted mass-ordering (IMO) of light neutrinos and thus different 
neutrino phenomenology. Hence, the present mechanism can be distinguished from the previous 
one [1] by the on-going and upcoming experiments on the neutrino oscillations [2] and neutrinoless 
double-beta decays [3]. 

The oscillation data from solar and atmospheric neutrinos, and from the terrestrial neutrino beams 
produced in the reactor and accelerator experiments, have measured two mass-squared differences 
(Amjji, Amli) and two large mixing angles (612, #23) to good accuracy [3] [5j . The two compelling 
features are [1] [5] : (i) the atmospheric neutrino mixing angle #23 has only small deviations from its 
maximal value of #23 — 45° ; (ii) the reactor neutrino mixing angle #13 is found to be small, having its 
allowed range still consistent with #13 = 0° at 90%C.L. Hence, the pattern of (#23, #13) = (45°, 0°) 
is strongly supported by the experimental data as a good zeroth order approximation. It is important 
to note that this pattern corresponds to the \i — r symmetry and Dirac CP conservation in the 
neutrino sector, where the \i— r symmetry is determined by both values of (#23) #13) = (45°, 0°) and 
the Dirac CP conservation is due to #13 = 0° . On the theory ground, it is natural and tempting 
to expect a common origin for all CP-violations, although the Dirac and Majorana CP-violations 
appear differently in the light neutrino mass-matrix of the low energy effective theory. Given such a 
common origin for two kinds of CP-violations, then they must vanish together in the [i— r symmetric 
limit. For the [l—t blind seesaw, we can uniquely formulate this common breaking in the Dirac mass 
matrix, leading to distinct neutrino phenomenology. 

With such a conceptually attractive and simple construction of the common breaking of two 
discrete symmetries, we can predict the fi—T breaking at low energies and derive quantitative cor- 
relations between the two small deviations, 23 — 45° and 13 — 0°, very different from that of the 
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previous NMO scheme pQ. Our predicted range of 8 13 can saturate its present experimental upper 
limit. The improved measurements of 823 will come from the Minos [7] and T2K experiments, etc, 
while #13 will be more accurately probed by the on-going reactor experiments, Daya Bay [10|[11| . 
Double Chooz[T2], and RENO [13], as well as the accelerator experiments T2K[8], NO^A [JJ] and 
LENA [15] . etc. We further derive the observed baryon asymmetry via leptogenesis at seesaw scale, 
and analyze the correlation between the leptogenesis and the low energy neutrino observables in 
the present IMO scheme. Especially, we deduce a lower bound on the reactor neutrino mixing an- 
gle $13 ^ 1° ; an d demonstrate that most of the predicted parameter space will be probed by the 
on-going Double Chooz, Daya Bay, and RENO reactor experiments. 

Finally, we will analyze the most general !L2®7t>2 symmetry structure of the light neutrino sector, 
and map it into the seesaw sector, where one of the Z2's is the fx— r symmetry Z^ T and another the 
hidden symmetry Z| (revealed in our recent work [T] for the NMO scheme) , which dictates the solar 
mixing angle 812 ■ We derive the physical consequences of the Z| for the most general light neutrino 
mass-matrix (without seesaw) and for the seesaw models (with different (jl—t breaking mechanisms). 
In particular, we analyze the partial violation of Z| in the presence of fi— r breaking for the ji— t 
blind seesaw, which leads to a modified new correlation between the \i — t breaking observables, very 
different from that of Ref. pQ. The determination of 812 is systematically studied for the current 
IMO scheme and the partial violation of Z| will be clarified. 

We organize this paper as follows. In Sec. 2 we present a unique construction for the common 
origin of the fj,— r and CP breakings in the neutrino seesaw with fj,— r blind right-handed neutrinos. 
Then, we give in Sec. 3 a model-independent reconstruction of light neutrino mass-matrix under 
inverted mass-ordering and with small fx — r and CP violations at low energies. In Sec. 4.1, we 
explicitly derive the low energy fj,— r and CP violation observables from the common breaking in the 
Dirac mass-matrix of the fx— r blind seesaw. These include the two small deviations for the mixing 
angles # 2 3 ~~ 45° and #13 — 0°, the Jarlskog invariant for CP- violations, and the M ee element for 
neutrinoless double-beta decays. In Sec. 4. 2 we study the cosmological CP violation via leptogenesis 
in our model, this can generate the observed baryon asymmetry of the universe. Using all the existing 
data from neutrino oscillations and the observed baryon asymmetry [161 [T7] , we derive the direct link 
between the cosmological CP-violation and the low energy Jarlskog invariant J . We further predict 
a lower bound on the reactor mixing angle 813 , and deduce a nonzero Jarlskog invariant J with 
negative range. We also establish a lower limit on the leptogenesis scale for producing the observed 
baryon asymmetry. In Sec. 5, we analyze the determination of solar mixing angle 812 and its relation 
to the hidden symmetry 7L\ in the light neutrino sector (without seesaw) and in the seesaw sector 
(with two different realizations of fx — r breaking). Finally, conclusions are summarized in the last 
section 6. 
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2. Common Origin of ft — r and CP Breaking from Neutrino Seesaw 
with Inverted Ordering 

The current global fit of neutrino data [1] for the three mixing angles and two mass-squared differences 
is summarized in Table- 1. We note a striking pattern of the mixing angles, where the atmospheric 
angle #23 has its central value slightly below the maximal mixing [18] of 45° and the reactor angle 
#13 slightly above 0° . So the neutrino data support two small deviations 623 — 45° and O23 — 0° of 
the same order, 

- 7.0° < (9 2 3 - 45°) < 5.5° , 0° ^ (#13 - 0°) < 9.5° , (2.1) 

at 90%C.L., with the best fitted values, (0 23 - 45°) = -2.2° and (0 13 - 0°) = 5.1°. This justifies 
a fairly good zeroth order approximation, 623 = 45° and 613 = 0°, under which two exact discrete 
symmetries emerge, i.e., the \x— r symmetry [19] and the Dirac CP conservation in the neutrino sector. 
It is clear that the \i — t symmetry and the associated Dirac CP-invariance are well supported by 
all neutrino data as a good zeroth order approximation, and have to appear in any viable theory 
for neutrino mass-generation. We also note that the #13 = 0° limit does not remove the possible 
low energy Major ana CP-phases, but since the Majorana CP-violation comes from a common origin 
with the Dirac CP-violation in our theory construction (cf. below), it has to vanish as the Dirac 
CP-violation goes to zero in the \i— r symmetric limit. 



Parameters 


Best Fit 


90%C.L. 


99%C.L. 


la Limits 


3<7 Limits 


Am^lO^eV 2 ) 


7.59 


7.26-7.92 


7.00-8.11 


7.39 - 7.79 


6.90-8.20 


Aml 1 (10- 3 eV 2 )(NM0) 


2.46 


2.26-2.66 


2.14-2.78 


2.34-2.58 


2.09-2.83 


Am 2 3 (lCT 3 eV 2 )(iM0) 


2.36 


2.18-2.54 


2.04-2.68 


2.25-2.47 


1.99-2.73 


#12 


34.5° 


32.8° - 36.0° 


32.1° - 37.2° 


33.5° -35.5° 


31.7° - 37.7° 


#23 


42.8° 


38.0° - 50.5° 


36.5° - 52.0° 


39.9° - 47.5° 


35.5° - 53.5° 


#13 


5.1° 


0°-9.5° 


0°- 11.3° 


1.8° - 8.1° 


0°- 12.0° 



Table 1: Updated global analysis [3] of solar, atmospheric, reactor and accelerator neutrino data 
for three-neutrino oscillations, where the AGSS09 solar fluxes and the modified Gallium capture 
cross-section [2D] are used. 

In our theory construction, we conjecture that all CP violations (both Dirac and Majorana types) 
have a common origin and thus they must share the common origin with the fi— r breaking. For the 
neutrino seesaw with heavy right-handed neutrinos blind to the fi— r symmetry, this common origin 
can only come from the Dirac mass-term. In the following, we first consider the minimal neutrino 
seesaw Lagrangian with exact fi— r and CP invariance, from which we will derive the seesaw mass- 
matrix for the light neutrinos. Diagonalizing this zeroth order mass-matrix we predict the inverted 
mass-ordering of light neutrinos and deduce the mixing angles, (823 , #13)0 = (45°, 0°) , as well as 
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a formula for the solar angle #12 ■ Then we will construct the common origin for the fi — T and CP 
breaking in the Dirac mass-matrix. Finally, we systematically expand the small fi—r and CP breaking 
effects in the seesaw mass-matrix to the first nontrivial order. 

2.1. fi — r and CP Symmetries of Neutrino Seesaw with Inverted Ordering 

The right-handed neutrinos are singlets under the standard model gauge group, and thus can be 
Major ana fields with large masses. This naturally realizes the seesaw mechanism [21] which provides 
the simplest explanation for the small masses of light neutrinos. For simplicity, we consider the La- 
grangian for the minimal neutrino seesaw [22j[23], with two right-handed singlet Majorana neutrinos 
besides the standard model (SM) particle content, 

£ ss = -LY t $£ R - L Yj& TV + -M^MrCM + h.c. 

= -l[M e £ R -17jrm D Af+ -M^MrCM + h.c. + (interactions) , (2.2) 

where L represents three left-handed neutrino-lepton weak doublets, I = (e, fi, t) t denotes charged 
leptons, u L = {u e , v , v T ) T is the light flavor neutrinos, and M = (N±, N2) T contains two heavy 
right-handed singlet neutrinos. The lepton Dirac-mass-matrix Mg = vYi/^/2 and the neutrino 
Dirac-mass-matrix m D = Y v arise from the Yukawa interactions after spontaneous electroweak 
symmetry breaking, ($) = (0, ) T 7^ , and the Majorana mass-term for Mr is a gauge-singlet. 
We can regard this minimal seesaw Lagrangian in Eq. (|2.2p as an effective theory of the general 
three- neutrino seesaw where the right-handed singlet N3 is much heavier than the other two (N\, N2) 
and thus can be integrated out at the mass-scales of (N\, N2), leading to Eq. (|2.2p . As a result, the 
minimal seesaw generically predicts a massless light neutrino [22] ; this is always a good approximation 
as long as one of the light neutrinos has a negligible mass in comparison with the other two (even if 
not exactly massless). Extension to the three-neutrino seesaw will be discussed in Sec. 4. 3. 

Let us integrate out the heavy neutrinos (N\, N2) in ()2.2[) and derive the seesaw formula for the 
3x3 symmetric Majorana mass-matrix of the light neutrinos, 

My ~ m D M R l m T D , (2.3) 

where m D is the 3x2 Dirac mass-matrix, and Mr is the 2x2 Majorana mass-matrix. The di- 
agonalization of M v is achieved by unitary rotation matrix U v via U^M V U V = D v with D v = 
diag(m 1 , m 2 , m 3 ) . 

The Lagrangian ()2.2[) is defined to respect both the fi— r and CP symmetries. Under the [i— r 
symmetry Zff", we have the transformation, Vn <-> pv T , where p = ± denotes the even/odd parity 
assignments of the light neutrinos under . Since the \i — t symmetry has been tested at low 
energy via mixing angles of light neutrinos, it is logically possible that the right-handed heavy 
Majorana neutrinos in the seesaw Lagrangian (|2.2p are singlets under (called "fi — t blind"), 
which is actually the simplest realization of fx— r symmetry in the neutrino seesaw. In this work we 
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consider that the right-handed Majorana neutrinos N to be fj,— r blind, i.e., both (Aq, N 2 ) are the 
singlets under Zg 7 ", and thus can be first rotated into their mass-eigenbasis without affecting the jjl— t 
symmetric structure of the Dirac mass-matrix m D . So, in the mass-eigenbasis of (N±, N 2 ), we have 
Mr = diag(Mi, M 2 ) ■ Under the [i — r and CP symmetries, the Dirac mass-matrix m D is real and 
obeys the invar iance equation, 



G T v m D 



D - 



with 



G„ 




(2.4) 



(2.5) 



Next, we note that due to the large mass-splitting of jjl and r leptons, the lepton sector can exhibit, 

in general, a different flavor symmetry Gi from the fi— r symmetry in the neutrino sector. The 

two symmetries Z 2 T and G^ could originate from spontaneous breaking of a larger flavor symmetry 

G f [2l] . Under the transformation of left-handed leptons F i £ G^ , we have the invariance equation 

rt 



of lepton mass-matrix, F}M»M\Fp 



MfMl . As we will show in Sec. 4. 2, we are free to choose an 



equivalent representation 



u}F e U e of 



from the start under which the left-handed leptons are 



in their mass-eigenbasis, where Ui is the transformation matrix diagonalizing the lepton mass-matrix, 



UjMgMjUg = Dj with D ( = diag(m e , m„, m T ) . This means that in the lepton mass-eigenbasis, 
the conventional Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix V [28] in the leptonic 
charged current (an analog of the CKM matrix [29] in the quark sector) is fixed by the transformation 
U v of neutrino mass-diagonalization, V = U u . We can further rotate the right-handed leptons into 
their mass-eigenbasis, without affecting the PMNS matrix, except making the lepton- mass-term 
diagonal in the seesaw Lagrangian (|2.2|) . i.e., Mi = diag(m e , m^, m T ) . 

Under the fJ,—T and CP symmetries, we find the Dirac mass-matrix m D to have the following 
form, 



m 



D 



(a 
b 



c 
c 



a 2 a 
a 2 c 
a 2 c ) 



(2.6) 



cj 1 a 
a 1 b 

\ M 

with all elements being real, and a 1 = y/ fh Q M\ , a 2 = \Jm^M 2 . As will be shown shortly, 
the parameter fh is defined at the seesaw scale and equals the nonzero mass-eigenvalue of the 
light neutrinos at zeroth-order under the fx — r symmetric limit. In (|2.6p we have also defined four 
dimensionless parameters, 

(a, b) , , , (a', c) 



(a, b) 



(a, c) 



V mo Mi ' \/ r moM 2 
Then, we find it convenient to define a dimensionless Dirac matrix, 

(a a' \ 



(2.7) 



m 



D 



m D (m M R ) 2 



b 



(21 



Substituting the above into the seesaw equation (|2.3p . we derive the fi — r and CP symmetric mass- 
matrix for light neutrinos, 

(a 2 + a' 2 



M„ 



1 rp 

m D M R m D 



m (m D m T D ) 



m 



V 



(2.9) 



ab + a'c ab + a'c\ 

+ c 2 b 2 + c 2 

b 2 + C 2 / 

which we call the zeroth order mass-matrix. In the next subsection we will further include the small 
[i— t and CP breaking effect. Note that from ()2.9|) . we have det(M u ) = 0, which generally holds in 
any minimal seesaw. 

Diagonalizing the mass-matrix ()2.9|) . we derive the mass-eigenvalues and mixing angles at zeroth 
order, 



m 



1,2 



mo 



mo 
2 

. 



(a 2 + a' 2 + 2b 2 + 2c 2 ) T Vi(a 2 + a' 2 ) - 2(& 2 + c 2 )] 2 + 8(ab + a'c) 2 



tan 2#i2 



2V2|a6 + a'c| 



723 



45 c 



713 



0°, 



(2.10a) 
(2.10b) 
(2.10c) 



| a 2 +a /2_ 2 (62 + c 2)| ' 

where we have made all mass-eigenvalues positive and the mixing angles (#12, #13, #23) within the 

range [0, l|] by properly defining the rotation matrix. (As shown in Table- 1, the solar angle 612 

is most precisely measured and its 3cr range is below 37.7°, so we always have 2#i2 < \ and 

tan 20i2 > .) The mixing angles (#23) #13) = (45°, 0°) are direct consequence of the [i—t symmetry, 

but this symmetry does not fix Q\% . Eqs. (|2.10b )-( |2.10b ) show that the mass-spectrum of light 

neutrinos falls into the "inverted mass-ordering" (IMO), rh<i > m\ S> 7713 . 

Table- 1 shows that the ratio of two mass-squared differences, ,^ m ^ 1 , <C 1 • Since for the minimal 

|z\m 31 1 

seesaw model with IMO, the equation det(M u ) = leads to fh 3 = , so the above ratio requires 
the approximate degeneracy m 1 — rh 2 to be a good zeroth order approximation as enforced by the 
neutrino oscillation data. So, we will realize the exact degeneracy fh 1 = rh 2 for the fi— r and CP 
symmetric mass-matrix (|2.9p . by imposing the relations for Eq. (|2.10aD . 



(a 2 + a' 2 ) - 2(b 2 + (?) = 0. 



ab + a'c = . 



(2.11) 



As will be shown in the next subsection, including the common origin of [i— r and CP breaking in the 
neutrino seesaw can produce small non-degeneracy between fh 1 and rh 2 at the next-to-leading order 
(NLO). Since the mass-parameter fh is introduced in (j2.7[) for defining the dimensionless parameters 
(a, b, c), we can now fix fh by defining 



m n 



mi 



??? 



2 ! 



(2.12) 



as the zeroth order mass-eigenvalue of light neutrinos, under the normalization condition, 

(a 2 + a' 2 ) + 2(6 2 + c 2 ) = 2. (2.13) 
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Combining this relation to Eq. (|2.1ip . we can deduce, 

a 2 = 2c 2 = 1 - 2b 2 , a' 2 = 2b 2 , c 2 = - - b 2 , a'c = -ab , (2.14) 

where we see that three of the four parameters, (a, a', c), can all be solved in terms of b . The last 
equation in (|2.14p is not independent, but it helps to fix a relative sign. We note that in (|2.9p the fi— t 
symmetric seesaw mass-matrix M u contains five parameters, the mass-parameter fh and the four 
dimensionless parameters (a, b, c, a'). The inverted mass-spectrum have imposed a LO condition 
m 1 = fh 2 , which results in two constraints in (|2.1ip . and the normalization condition m = rh 1 in 
(|2,12p leads to the third constraint (|2.13p . In consequence, we end up with only two independent 
parameters, rh and b . 

We note that under the condition of (|2.1ip , the mixing angle 812 given by (|2.10cp has no definition 
at the zeroth order (the fi—T symmetric limit) due to the vanishing numerator and denominator in 
the formula of tan 26 12 ■ But including the small fi — T breaking effect will generate the nonzero 
expression of B\i at the NLO even though its final formula does not depend on the fi — r breaking 
parameter (cf. Sec. 2. 2). As we will show in Sec. 2.2, the fi — r breaking arises from deviation in 
the element c of m D , so we can apply the 1 'Hopital rule to the expression of tan 28±2 by taking the 
first-order derivatives on its numerator /denominator respect to c and deduce, 

un2e ^^k = ^k' (2 ' 15> 

which is consistent with (|4,5p of Sec. 4. 1 from the explicit NLO analysis. For the case with fi — T 
breaking arising from deviation in the element b of m D , we can apply the 1 'Hopital rule again to 
infer the formula, 

tan20 12 = J4- = ^M, (2.16) 
V2\b\ \a\ 

which is the inverse of \2.15\) . As will be shown in Sec. 5. 2, the different forms of jjL— r breaking will 
affect the determination of the solar mixing angle 0\2- But it is worth to note that the expression 
of #12 is fixed by the jjl— t symmetric m D as in (j2. 15[) or (|2.16p . and does not explicitly depend on 
the fi — T breaking parameter. We will systematically analyze these features in Sec. 5 and clarify the 
difference from our previous construction [1] . 

2.2. Common Origin of fi — t and CP Breaking in the [i — T Blind Seesaw 

In this subsection, we will construct a unique breaking term providing a common origin for both fi—T 
and CP breaking. From this we will further derive predictions of the common fi—T and CP breaking 
for the low energy light neutrino mass-matrix, by treating the small breaking as perturbation up to 
the first nontrivial order (Sec. 4). We will analyze the seesaw-scale leptogenesis and its correlations 
with the low energy observables in Sec. 4. 2. 

As we have explained, the fi — T symmetry serves as a good zeroth order flavor symmetry of 
the neutrino sector, which predicts #13 = and thus the Dirac CP-conservation. Hence, the fi — T 
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symmetry breaking is generically small, and must generate all Dirac CP-violations at the same time. 
On the theory ground, it is natural and tempting to expect a common origin for all CP-violations, 
even though the Dirac and Major ana CP-violations appear differently in the light neutrino mass- 
matrix of the low energy effective theory. For the two kinds of CP-violations arising from a common 
origin, then they must vanish together in the fi— r symmetric limit. 

Different from our previous study pQ, we consider the heavy right-handed neutrinos to be jjl— t 
blind in the neutrino seesaw. Thus the Majorana mass-matrix M R of the right-handed neutrinos 
must be li—t singlet. Hence, we deduce that the unique common origin of the li—t and CP breaking 
must arise from the Dirac mass-matrix of the seesaw Lagrangian \2.2\) . For the minimal seesaw, the 
most general form of m D is 



(2.17) 



where the scaling factors a 1 = y^m^Mi and o 2 = y / fh Q M2 are real mass-parameters as defined in 
Eq. (|2.6p . The six elements of m D can be complex in general. But there are three rephasing degrees 
of freedom for the left-handed lepton-doublets. So we can always rotate the three elements in the 
first column of m D to be all real, hence the remaining CP phases (associated with the li—t breaking) 
have to appear in the elements c 1 and c 2 because a' cannot break \i—t symmetry and thus should 
be real. We have conjectured that all CP violations arise from a common origin, which then must 
originate from the li — t breaking; so we can formulate such a common origin as a single phase in 
either c 1 or c 2 in the minimal construction, where the other two elements in the second column 
of m D should be real. Hence, we present a unique minimal construction to formulate the common 
origin of Li— r and CP breaking in the Dirac mass-matrix m D as follows, 

/ cJi a a 2 a' \ 





( 


a 


a' 


\ 




( 


°1 


a 


cr 2 a 1 


\ 


m D = 




\ 










°l 


bi 


a 2 c x 






\ 




~ J 2 


J 




\ 


°l 


b 2 


<?2 c 2 


/ 



m 



D 



2.18) 



a x b cr 2 c(l-C') 
\a x b a 2 c(l-(e^)J 

where the dimensionless parameters — 1 < £' < 1, 0^C<1, and the CP-phase angle uj 6 [0, 2ir) . 

b since (& l5 b 2 ) are already made real and thus cannot serves as the 



Here we have set b 



common source of the fx— r and CP breaking. Inspecting (|2. 18j) we see that, for any nonzero £ and 
uj, the li — t and CP symmetries are broken by the common source of Qe lul . We could also absorb 
the real parameter £' into c by defining d = c(l — CO • Thus we have, 

\ 



<7i a 



( » 

i 



a x b 



o~ a 



CJn c 



(2.19) 



with 



C'e 1 ' 



a 2 c'{l-C"e^) ) 
(e iuJ - (' 

i-C ' 
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(2.20) 



Given the ranges of (£, £') as denned above, we see that the corresponding new parameter C," of the 
fx — t breaking has a much larger range, including values within 1 < < 3 (when |C| 5 IC'I ^ 0.6 
for instance), which are beyond the perturbative expansion. We find that if enforce < 1 , the 
parameter-space of (j2. 19[) becomes smaller than f)2 . 18|) and insufficient for making the model fully 
viable. This means that our formulation of (|2.18p is more general and has larger parameter-space for 
making theoretical predictions. Hence, we will apply (|2.18p for the physical analyses below. 
We note another formulation of such a breaking in the Dirac mass-matrix m D , 

/ cr 1 a a 2 a' \ 
m D = a x b <r 2 c(l - (e lu >) , (2-21) 

[a.b a 2 c(l-C0 / 
which is connected to (|2.18p by a f/,— r transformation for the light neutrinos v = (u e , v , v T ) into 
v' = (z/ e , v T , i/^) , via v = G v u' , with G u [p = 1] defined in Eq. (|2.5p . Accordingly, the mass-matrix 
()2.2ip transforms as, 

m D — > m' D = G^m D = m D , (2.22) 

which goes back to (|2.18|) . So the two different formulations (|2.18|) and (|2.21|) just cause the [/,— r 
asymmetric parts in the seesaw mass-matrix M v = m D M^ to differ by an overall minus sign. As 
we will comment further in Sec. 4. 1, this does not affect our predictions for the physical observables 
and their correlations. So we only need to focus on the formulation (|2.18p for the rest of our analysis. 

We may also first rotate the three elements in the second column of (|2.17|) to be real and then 
formulate the common origin of /U — r and CP breaking as follows, 



m 



D 



{ <J 1 a 

M(W) 

V ^(l-Ce*") 



<T 9 C 



(To C 



(2.23) 



/ 



As will be clarified in Sec. 5, this will lead to the determination of solar mixing angle #12 as in 
(|2,16p . in contrast to (|2.18p which predicts a different 612 as in ()2.15|) . Here 612 is explicitly fixed 
by the fx— r and CP symmetric parameters of m D in either case. But, we find the predictions for all 
other [i— t and CP breaking observables and their correlations to remain the same as those from the 
construction in (|2. 18j) . 

Finally, it is interesting to note that for an extended Higgs sector (consisting of two Higgs dou- 
blets or more) we can generate all CP-phases in the Dirac mass-matrix m D via spontaneous CP 
violation [30], which is beyond the current scope and will be elaborated elsewhere [31j . 



2.3. Perturbative Expansion for \x — r and CP Breaking 

Let us first consider the 3x3 mass-matrix M v light neutrinos, which can be generally presented as, 

A B\ B 2 \ / A B B \ ( 5 A 8B1 5B 2 



My 



Ci D 

C 2 



\ 



Do 
Co 



+ 



V 



5d 5D 

sc 2 



11 



M^ + 5M U = M^+5MP+0(Cf) 



(2.24) 



where the zeroth order matrix My corresponds to vanishing fi — T breaking with Q = , and the 
NLO mass-matrix 5M^ includes the \i — r breaking to the first nontrivial order. We find it useful 
to further decompose SM^ into the fi — t symmetric and anti-symmetric parts, 



5MP = 8M' + 6MS = 



( 5 A 5B S SB, 
SC S 5D 

so, 



/ 6B a -5B a 

sc a 

V -5C a 



(2.25) 



with 



SB, 



1 



{5B l + 5B 2 



SB a 



5C S = -(5C 1 + 5C 2 ) 



1 



1 



(SB X -SB 2 ) , 



SC a = - (6d -SC 2 ) . 



(2.26a) 
(2.26b) 



This decomposition is actually unique. 

Prom our construction in the previous subsection, the // — r and CP breaking Dirac mass-matrix 
m D as well as the Majorana mass-matrix M R is uniquely parameterized as follows, 

f a l a a 2 a' \ 

M R = diag(Mi, M 2 ), 



m 



D 



a x b a 2 c x 
\a x b a 2 c 2 ) 



(2.27) 



with o x 2 = \/^oMip and 

c x = c(l-C')» c 2 = c(l-Ce^). 

Thus, we can explicitly derive the seesaw mass-matrix for light neutrinos, 

/a 2 + a' 2 ab + a'c x ab + a'&A 



(2.28) 



M„ 



m 



V 



(2.29) 



¥ + cf & 2 + c x c 2 
b 2 + c 2 J 

Since the neutrino data require the \x — r breaking to be small, we can further expand M v in terms 
of small breaking parameter £ as, 



with 



M,(°) 



M u = 


mV* + 6M v = M(°) + SMp + 0(C 2 ) , 




(2.30) 




f a 2 + a' 2 ab + a'c afr + a'c^ 




n o 






= ™o 


& 2 + c 2 b 2 + c 2 


= ™0 


i 

2 


!)■ 


(2.31a) 




^ fr 2 + C 2 y 






2 / 






/ -a'c C' -a'c (e iuJ \ 










= "*o 


-2 c 2 C' -c 2 (C' + Ce ia; ) 








(2.31b) 




^ -2c 2 (e iu > j 
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where we have used the solution ()2.14p for the second step of (|2.31ap and the //— r breaking expression 
()4.58bj) for deriving ()2.31b[) . For our current model with the expansion up to 0{Q, £'), we deduce 
from P^Taj) - (^3Tb|l and (p^5l) - (p^6l) . 

A = m (a 2 + a' 2 ) = m , (2.32a) 
= m (a6 + ac) = 0, (2.32b) 
Co = D = m (k 2 + c 2 ) = -fh Q , (2.32c) 



and 



£A = 0, 5D = -m c z (C + (e tuJ ) , 

5B S = -±m a'c((' + (e l "), 5C S = -m c 2 (C' + (e M ) , ( 2 . 33 ) 
<yB a = -lm a'c(£ - Ce w ) , 5C a = -m c 2 ((' - Qe^) . 



2 

Note that from (|2.33|) we can compute the ratio 



SB a a' b 



(2.34) 



5C a 2c a ' 

where in the last step we have used the resolution (|2.14|) . It is interesting to note that the ratio 
(|2.34p of the fi— r asymmetric parts in the light neutrino mass-matrix M v only depends on the [i — r 
symmetric elements of the Dirac mass-matrix m D . This ratio just corresponds to the determination 
of the solar angle #12 in (|2.15|) and will be further confirmed later by the full NLO analysis of Sec. 4. 1. 



3. Inverted Ordering: Reconstructing Light Neutrino Mass Matrix 
with /x — r and CP Violations at Low Energy 

In this section, we give the model-independent reconstruction of the Majorana mass-matrix for 
light neutrinos under inverted mass-ordering (IMO), in terms of the low energy observables (mass- 
eigenvalues, mixings angles and CP phases). We expand this reconstruction by experimentally well- 
justified small parameters up to the next-to-leading order (NLO). Applying this reconstruction formu- 
lation to our model will allow us to systematically derive the physical predictions for the correlations 
among the low energy observables as well as for the link to the baryon asymmetry via leptogensis at 
the seesaw scale. 



3.1. Notation Setup and Model-Independent Reconstruction 

Let us consider the general 3x3 symmetric and complex Majorana mass-matrix for the light neutrinos, 

/m ee m efl m er \ I A B x B 2 \ 

M v = m w = C\ D (3.1) 

m TT J 



\ 



\ 



C 2 J 
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In the mass-eigenbasis of charged leptons, the neutrino mass-matrix M v can be diagonalized by a 
unitary transformation V(= U v ) , i.e., V T M U V = D v = diag(m 1 , m 2 , m 3 ) , and thus we can write 
the reconstruction equation, 



M v = V*D V V* . 



(3.2) 



The mixing matrix V can be generally expressed as a product of three unitary matrices including a 
CKM-type mixing matrix U plus two diagonal rephasing matrices U' and U" , 



V = U"UU' , 



( 



U 



" SaCnr. 



Ss^a C s S a S x € D C- s C a -\- S s S a S x C 



-iS D 



S x 6 



\Ss$a H~ C s C a S x G D C s S a S s C a S x C D C a C x 

U' = diag(e^i, e^, e^s), u" = diag(e iQ i, e ia *, e i 



(3.3a) 
(3.3b) 

(3.3c) 



where 5 a is the Dirac CP-phase. For notational convenience, we have denoted the three neutrino mix- 
ing angles of the PMNS matrix as, (6 12 , 6 23 , 13 ) = {6 S , a , X ) , by following Ref. [23]. We will fur- 
ther use the notations, (s s , s a , s x ) = (sin# s , sin6> a , sin^.) and (c s , c a , c x ) = (cos9 s , cos0 a , cos9 x ) . 
For the diagonal rephasing matrix U' , only two of its three Majorana phases are measurable (such as 
4>3 — (pi and <f>2 — 4>i) after extracting an overall phase factor. The matrix U" contains another three 
phases which associate with the flavor-eigenbasis of light neutrinos and are needed for the consistency 
of diagonalizing a given mass-matrix M v . 

For convenience we define the rephased mass-eigenvalues D u = U'*D U U'^ = (mi, TO2, fh 3 ) = 
(m ie - i2 *i, m 2 e- i2 ^, 

m 3 e_ * 2< ^ 3 ) ' so the reconstruction equation ()3.2|) becomes, 



M v = V'*D V V' ] , ( V' = U"U ) 



(3.4) 



Thus, we can fully reconstruct all elements of M u in terms of the rephased mass-eigenvalues (mi, fhs), 
the mixing angles (8 S , 6 a , 9 X ), the Dirac phase So, and the rephasing phases ctj (which do not appear 
in physical PMNS mixing matrix), 



m. 



??? 



/'/' 



m^ 



e -i2 Ql 



e -i2a 2 



2 2' — : I 22' — I 2 — 2i<5^> 
CgC^Tfi^ -\- S s C x Tfi<2 ~\~ S x & 



2 2, 



■nu 



(s s c a - c s s a s x e z5D ) 2 m 1 + (c s c a + s s s a s x e l5D ) 2 m 2 + s 2 c 2 m 3 
(s s s a + c s c a s x e tbu ) 2 rh x + (c s s a - s s c a s x e l5D ) 2 fh 2 + c 2 c 2 m 3 



(3.5a) 
(3.5b) 
(3.5c) 



777 



e/i 



m„ 



77? 



-i( ai +a 2 ) 



-*(«i+«3) 



-i(a 2 +a 3 ) 



CsC x (s s c a c s s a s x e d )tti^ — s s c x (c s c a -\- s s s a s x e D )m 2 -\-s a s x c x e D ni^ , (3.5d) 
CsC x (s s s a +c s c a s I ,e D jiTJj — s s c x (c s s a — s s c a s x e D )m 2 — c a s x c x e D m 3 , (3.5e) 

(^sQj CgSaS^e D )(^sS a "I - C s C a S x e D )?7lj 
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+ {c s c a + s s s a s x e i ^ D )(c s s a — s s c a s x e^ D )m2 — s a c a c^.m^ , (3.5f) 

where among the Majorana phases </>i 2 3 (hidden in the mass-parameters m l 2 3) only two are 
independent because an overall phase factor of U' can be absorbed into the diagonal rephasing- 
matrix U" . For the case with a vanishing mass-eigenvalue (such as m 3 = in our present model), 
only one independent phase combination, say e*^ 2 "^ 1 ^, will survive. If we impose \i—t symmetry on 
the light neutrino mass- matrix M u , we can deduce [I], 



'at v xJq 



(45°, 0°), « 20 = «3o- (3-6) 



The solar mixing angle 9 S is independent of the [i — r symmetry and is thus left undetermined. To 
predict 6 S , we will uncover a new flavor symmetry beyond the Z 2 r (cf. Sec. 5). 

3.2. Reconstruction of Light Neutrino Mass Matrix with Inverted Ordering 

Now we are ready to apply the above general reconstruction formalism to the inverted mass-ordering 
(IMO), m 2 > m 1 S> m 3 , with m 3 = (as predicted by the present minimal seesaw model), in 
contrast to our previous model which predicts the normal mass-ordering (NMO) pQ . We introduce a 
small mass-ratio for light neutrinos, 



m 2 ~ m i Am 2 
m? Am? 



y = 1 2 1 = -j— f- = 0.029 - 0.036 < 1 , (3.7) 

"^13 



as constrained by the neutrino data at 90%C.L. (Table-1). So it is sufficient to make perturbative 
expansion in y' up to its linear order. Thus, at the zeroth order of y' , we have equal mass-eigenvalues, 
m 10 = m 20 = m . Under the y'-expansion up to next-to- leading order (NLO), m i = m + 5m i , we 
have 

, ^ 2(5m 2 - dm^ = 2(m 2 - 

We can define another small ratio z = - — - = 0(y') , and deduce, 

m 1 

5m 1 = zmi, 5m 2 = ( z + — J m 1 , (3-9) 



where m l = y Am^ is fixed by the neutrino data, and m = m 1 — 5m 1 = (1 — z)m l ~ v^Ar 

Next, we consider the mixing angles and CP-phases. Since the neutrino oscillation data strongly 
support the /i— r symmetry as a good approximate symmetry (|3.6|) , we can define the small deviations 
from the general fi— r symmetric solution (|3.6p . 

s a = Q a -j, s x = e x -o, (3.10) 

which characterize the fjr-r symmetry breaking. From the data in Table-1, we can infer the constrained 
% C.L. ranges, 

< Si ^ 0.027, < 6l < 0.015. (3.11) 
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For our analysis we will systematically expand the small parameters (<5 a , S x , y' , z) up to their linear 
order. For the Majorana CP-phases, (j) 3 drops due to m 3 = ; we also remove an overall redundant 
Majorana phase <j) l (from U') into the redefinition of a- (in U"). So, the remaining independent 
Majorana phase is only (j> , 



a j = a i + ^i' (i = 1 ; 2 ; 3 ) 



+ 



(3.12a) 
(3.12b) 



The expansion up to the NLO for our current reconstruction analysis will include {5a 1: 5a 2 , 5a 3 , 5<t>) ■ 
The solar angle 6 S (= #12) is independent of the fM—T breaking and thus receives no NLO correction. 
Furthermore, we note that the Dirac phase e^ D is always associated with the small mixing parameter 
Sx (— $x) ) so h only appears at the NLO and thus receive no more correction at this order of 
expansion. 

Finally, we give a summary of all relevant NLO parameters in our reconstruction analysis, 



(y', z , 5a, 6 X , <5«2> 6a 3 , 



(3.13) 



Each of them is defined as the difference between its full value and zeroth-order value under the [l—t 
symmetric limit. In Sec. 4 we will derive these deviations from our seesaw model for the common 
origin of fi— r and CP breaking, and analyze their correlations. 

Making the perturbative expansion of f)3. 13|) under the inverted mass-ordering, we first deduce 
the LO form of the light neutrino mass-matrix (|3.ip . 



m (o) =A = m e- 2iai ° ( cl + a^e"* 2 * 



1 

1 



-i(aio + 



Q20) f 



1 - e 



-i2<t>o 



Dn 



(3.14a) 
(3.14b) 

(3.14c) 



where we have also matched to our notation of in flZjUP . Then, we derive elements of the 

NLO mass-matrix SMj, from 



dm,™ 



"'"'fiH 



5m& 



SA = m e- im ™ 
SB t 



8mi l ) = 5Bo 



z + -±yf - i2(s 2 s 6(f) + Sax) 



^ e -*(<5l0+°20) 



V2 

!!!£ e -*(°io+«2o) 



V -e i5D S x + i2c s s s 



^y> + e iS °S x + i2c s s s 



SC\ = rriQe 



SC2 = moe 



-i2a 



-i2a 



I + -fy' - 6 a - i{c 2 s S<j> + Sa 2 ) 



z c; 



- + jy' + 6 a - i{cj 54> + 5a 3 ) 



(3.15a) 
(3.15b) 
(3.15c) 

(3.15d) 

(3.15e) 
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8m$ = 



8D 



| + jl/ ~ 8<t> + 8a 2 + 8a 3 ) 



where we have matched to our notation of 8M„ as defined in (|2.24|) , In the above formulas, we 



(3.15f) 



'aO, 



Ixo) = (f , 0) and a 



20 



a 



30 ' 



have used the fi — T symmetric relations for the LO parameters 
as well as 711,3 = . 

From (|2.25p . we can uniquely decompose the elements of 8M„ in (|3. 15|) as the fi — r symmetric 



and anti-symmetric parts, 8My 



8M^ + 8M£ , with their elements given by, 



5B S = 

5B a = 

8C S = 

8C a = 



8B X 


+ 8B 2 




2 


8B X 


-8B 2 




2 


8d 


+ 8C 2 




2 


Sd 


-8C 2 



mo 
7=2 



-i(a 10 +a 20 ) 



Cq S c. 



y + i2c s s s 



!^0 e -i(a 10 +a 20 ) e iSn § 
V2 



m e i2<5 2° 



+ 802 + 8a>3 



(3.16) 



8D. 



moe 



8a + ^(80:2 



8 a?. 



With these, we will be ready to apply the above reconstruction formulas (|3,14p . (|3.15p and (|3. 16|) to 
match with (|2.24p in our seesaw model at the LO and NLO, respectively. We will systematically solve 
these matching conditions in the next section, which allows us to connect the seesaw parameters to 
the low energy neutrino observables and deduce our theoretical predictions. 

For matching the seesaw predictions to our reconstruction formalism, we note that the latter was 
presented at the low energy scale so far. We need to connect the low energy neutrino parameters to 
the model predictions at the seesaw scale, where the possible renormalization group (RG) running 
effects should be taken into account in principle. Such RG effects were extensively discussed in 
the literature [25J , and can be straightforwardly applied to the present analysis. Below the seesaw 
scale, heavy right-handed neutrinos can be integrated out from the effective theory and the seesaw 
mass-eigenvalues m-{j = 1,2,3) for light neutrinos obey the approximate one-loop RG equation 
(RGE) [25], 

dm- 



dt 



a 



(3.17) 



to good accuracy [26] , where t = ln(/i//i ) with \i the renormalization scale. For the SM, the 
coupling-parameter S~ — Zg\ + 6?/f + A , with (g 2 , y t , A) denoting the SU(2)l weak gauge coupling, 
the top Yukawa coupling and Higgs self-coupling, respectively. Hence, we can deduce the running 
mass-parameter m- from scale /i to fi , 



exp 



1 



16tH 



a(t') dt' 



mAii ) 



(3.18) 



with t = ln(/j/ / u ) . In the present analysis we will choose, (/x , fi) = (Mz, Mi) , with Z boson mass 
Mz representing the weak scale and the heavy neutrino-mass M\ characterizing the seesaw scale. 
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Consider the minimal neutrino seesaw with inverted mass-spectrum, m 2 •> m 1 3> m 3 = . We 
note that the zero-eigenvalue m 3 and the mass ratio y' do not depend on the RG running scale \x . 
So we can derive the running of the two nonzero mass-parameters from weak scale to seesaw scale, 

m x = m 1 (M 1 ) = Xim^Mz), (3.19a) 
m 2 = m 2 (Mi) = Xim 2 (M z ) = \A + y' m x , (3.19b) 

with Xi = x(Mi,Mz)- In Sec. 4, we will compute the RG running factor X\ = x(MljMz) 
numerically, which depends on the inputs of initial values for a 2 = g 2 /(A7r) , y t and the Higgs boson 
mass Mh , via the combination a defined above. Using the electroweak precision data [171 127] . 
a ~ \M Z ) = 29.57 ± 0.02, m t = 173.1 ± 1.4 GeV, and the Higgs-mass range 115 < M H < 149 GeV 
[90% C.L.] for the SM, we find the running factor X {Mi,M z ) ~ 1.3 - 1.4 for Mi = 10 13 - 10 16 GeV. 
Other running effects due to the leptonic mixing angles and CP-phases are all negligible for the present 
study since their RGEs contain only flavor-dependent terms and are all suppressed by y^ = 0(10 ) 
at least [25]. For the analyses below (Sec. 4), we will first evolve the mass-parameters from the seesaw 
scale Mi down to the low energy scale for neutrino oscillations, and then match them with those 
in our reconstruction formalism. Including such RG effects just requires to replace the light mass- 
eigenvalues (m l5 m 2 ) at seesaw scale Mi by the corresponding (m 1 , m 2 ) at low energy, and vice 
versa. 

4. Predictions of Common fjt — r and CP Breaking with Inverted Ordering 

In this section we apply the reconstruction formalism (including the RG running effects) in Sec. 3. 2 to 
our common fx— t and CP breaking seesaw in Sec. 2. 3. Then, we systematically derive the predictions 
for the low energy neutrino observables. This includes the nontrivial correlation between two small 
fi—r breaking parameters 8 X (= $13 — 0) and 5 a (= O23 — \ ) ■ Furthermore, we study the correlations 
of #23 — 45° and #13 with Jarlskog invariant J and neutrinoless /3/3-decay observable M ee . Finally, 
we study the matter-antimatter asymmetry (baryon asymmetry) via leptogenesis in the fJ,—T blind 
seesaw, and establish the direct link with low energy neutrino observables. Furthermore, we will 
derive a nontrivial lower bound on the reactor mixing angle, #13 > 1° , and restrict the Jarlskog 
invariant into a negative range, —0.037 < J < —0.0035 . 

4.1. Predicting Correlations of Low Energy Neutrino Observables 

Both fj, — t and CP violations arise from a common origin in the seesaw Lagrangian of our model, 
which is characterized by the breaking parameter ^e luJ and shows up at the NLO of our pertur- 
bative expansion. Hence, in the light neutrino mass-matrix, the small // — r breaking parameters 
(^a; S x ) together with all CP-phases are controlled by £ and uj . In the following, we will use the 
reconstruction formalism (Sec. 3. 2) under IMO for diagonalizing the light neutrino mass-matrix at 
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the NLO. Then, we will further derive quantitative predictions for these low energy observables and 
their correlations. 

We first inspect the reconstructed LO mass-matrix M„ in (|3.14p . Matching (|3. 14j) with our 
model prediction (|2.31ap at the same order, we find the solutions, 



a 



10 



a 



20 



m 



in 



??? 



20 



Po = 0, 
mo, m 3 = , 



a 2 = 2c 2 = 1 - 2b 2 



2b 2 



1,2 I 

b , a c 



-ab , 



(4.1a) 
(4.1b) 
(4.1c) 



which is also consistent with Eq. (|2.14|) . Here all the LO CP-phases (a 10 , a 



10' "20' vo) 



because the 



original CP-violation in the seesaw Lagrangian vanishes in the Q = limit (Sec. 2. 2). 

Then, we analyze the NLO light neutrino mass-matrix 5Mu , as given by (|2.25p of our model 
and by the reconstruction formula (|3.15p . We match the two sets of equations at the low energy for 
the \x — t symmetric elements, 



6A = 



5B„ 



■m, 



°a'c((' + (e lu ) 



V2 



z + ^-y' -i2(s 2 5(j) + 8ai) 



Ca S s 



y' + i2c s s s 



z + y?/ - i(2c s 5(f) + 5a 2 + 5a 3 ) 



5D, 



and for fi — r anti-symmetric elements, 



8B a -- 

sc a 



in 



°a'c(('-(e l 



m Q c 2 ((>-(en 



V2 



5 a + -^{5a 2 - 5a 3 ) 



(4.2a) 
(4.2b) 

(4.2c) 

(4.3a) 
(4.3b) 



where using Eq. (|3.19p we have run the mass-parameter fh Q from the seesaw scale down to the 
corresponding m at low energy for the left-hand-sides of Eqs. (|4,2p and (|4.3p . 

From the [i— r symmetric Eqs. (|4.2b )- (|4.2b ). we can infer six independent conditions for the real 
and imaginary parts of (5 A, 5B S , 5C S ) , respectively, 



5a 1 



CsS s i 

y 



V2 
2V2 c„a 



2 



1 



Z C, , 

- + — y 

2 4 y 

(2c 2 54> + 5a 2 + 5a 3 ) 



2 

■2 



i'c (C' + C cosoj) , 
-a'cCsina; , 



- 2 (C + c 

-c 2 C sin u , 



cosw 



(4.4a) 
(4.4b) 
(4.4c) 

(4.4d) 
(4.4e) 
(4.4f) 
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Thus, with the aid of (|4.4ap we take the ratio of (|4.4c|> and (|4.4ep . and derive 

a' V2b , . 

tan26» s = =- = , (4.5) 

V2c a 

which coincides with (|2.15p in Sec. 2. 1. Using Eq. (|4.5|) . we deduce from Eq. (|4.1cp . 

a = p a cos29 s , b = p a —=sm2d s , (4.6a) 

v2 

a' = p a , sin 26 s , c = —p a , —=■ cos 29 s , (4.6b) 

v2 

with p , p , = ± denoting the signs of (a, a'). Here we see that the four dimensionless LO param- 
eters (a, a', b, c) in the Dirac mass-matrix (|2.18p are fixed by the solar mixing angles 6 S , since the 
conditions in ()4.1cp make three of them non-independent. Finally, we further resolve (|4.4p and derive 
the NLO parameters, 

y = -2cos26> s ((' + (cosu) , (4.7a) 
z = s 2 s cos26 s (C' + Ccosw) , (4.7b) 
6^ = -\s 2 s {c 2 s - s 2 s )( sinw, (4.7c) 

5<f> = ^(c 2 - S 2 )( S mu, (4.7d) 

8a 2 + Sa 3 = s s (s 2 — c s ) £ sin u; . (4-7e) 

It is interesting to note that the present model predicts a generically small Major ana CP-phase 
angle at low energy, <fi = 5(f) = 0(Q , in contrast to our soft breaking model [T] where the low energy 
Major ana CP-phase angle (^23) is n °t suppressed. 

Next, we analyze the fi—r anti-symmetric equations (|4.3ap - (j4.3bp for SMy . With (|4.6[) . we can 
deduce from (|43a>(|43E) . 

-sin20 s cos20 s (C' - (e iu ) = -e i&D 8 x , (4.8a) 



X - cos 2 26 8 (C' - Ce^) = 5 a + ^(<5a 2 -<fe 3 ), (4.8b) 



which decompose into 



cos^d^x = — - sin 2# s cos 20 s (£' — £ cos u;) , (4.9a) 

sin^D $ x = - sin 20 s cos 20 s (£ sin u;) , (4.9b) 

J a = icos 2 2^(C'-Ccosw), (4.9c) 

<fe 2 -5a 3 = -cos 2 26» s (C sin a;) . (4.9d) 

Thus the Dirac CP-phase angle 5d can be derived from the ratio of (|4.9ap and (|4.9blh 

C sino; 5a 2 - Sa 3 

tan<5 D = = — . (4.10) 

Q COS U) — Q 2d a 
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With Eqs. (pTTafl . (KlQ\\ and (021, we finally deduce, 



, / . I i , . 2 tan /,-,,% 



and thus 



sin 20 sin 20 

cos5 D 5 x = ^ (y' + 4cos20 s C') = — ^ {y' + 4 cos 20 s cos u C) , (4.12a) 

S a = ^k( y ' + 4 cos 20 s C) = -^i(y' + 4cos20 s co S a;C), (4.12b) 

<5a 2 — 5a 3 = 2tan<5£><5 a . (4.12c) 

Prom Eqs. (|4.12ap and (|4.12bp . we derive a nontrivial correlation between the low energy fj, — r 
breaking observables 5 a and 5 X , 

5 a = — cot 26 s cos 8d 5 X . (4-13) 

This shows that at the NLO the two small \i— t breaking parameters are proportional to each other, 
5 X oc 5 a . Because of | cos 5d\ ^ 1 , we can infer from Eq. (|4,13p a generic lower bound on 5 X , for 
any nonzero 5 a , 

6 X ^ |5 a |tan20 s , (4.14) 

where we have 5 X = Oyj G [0, f ] in our convention. It is worth to note that our previous soft breaking 
model [T] also predicted a correlation and a lower bound, 

S a = — cot 6 S cos <5d 5 x , (Prediction of Ref. [1]), (4.15a) 
=> 5 X ^ |5 a |tan0 s , (4.15b) 

where the quantitative difference from the present predictions is that we have the coefficient cot 29 s 
in Eq. (|4. 13j) as compared to cot0 s in Eq. (|4.15a|) . In fact, this is a profound difference. From the 
present oscillation data in Table-1, we observe that the deviation of the solar angle 9 S (= 612) from 
its maximal mixing value is relatively small, 

9.0° < 45° -0i2 < 12.2°, (at90%C.L.), (4.16) 

and this limit only relaxes slightly at 99% C.L., 7.8° < 45° — 0i2 < 12.9° . Hence, we see that the 
range of the deviation 45° — 0i2 is at the same level as the two other small deviations 023 — 45° and 
#13 — 0° shown in Eq. (|2.ip . So, we can define a new naturally small quantity, 

S s = j -0s, (4.17) 

and make expansion for 5 S as well. Then, we immediately observe a qualitative difference between 
cot 20 s ~ 25 s <S 1 in (|4.13p and cot S ~ 1 + 25 s > 1 in (|4.15aD . Hence, we can rewrite the two 
correlations (|4.13|) and (|4.15a|) in the well expanded form, 

5 a ~ — 2 cos Sd (S s 5 x ) ^ 5 X , (Current Prediction), (4.18a) 
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5 a ~ - cos S D 5 X = 0{8 X ), (Prediction of Ref. pQ). (4.18b) 



Two comments are in order. First, we deduce from (|4.18p the following patterns of the three mixing 
angles, 

(TT TT \ 

— - S s , — - 0(5 S 5 X ), 5 x j , (in the current model), (4.19a) 

(012, #23, 0is) = (| - 5 s ,j - 5 a , 4) , (in the model of Ref. P), (4.19b) 

where for the current model Eq. (|4.19ap predicts a nearly maximal atmospheric angle #23 — \ ; while 
for the soft-breaking model [T], Eq. (|4. 19b|) allows all three deviations to be comparable. Second, for 
each given nonzero 5 a = 623 — \ , we can deduce the lower limits on S x = 613 from (|4.18p . 

I S I 

5 X ^ -L^- > \6 a \, (Current Prediction), (4.20a) 
25 s 

5 X ^ \S a \, (Prediction of Ref. [I]). (4.20b) 

Given the 99% C.L. range of 7.8° < S s < 12.2° , we derive the lower limit from (gUP or (|4.20ap for 
the present model, 

013 > (3.6 ~ 2.1)|023 - 45° I, (4.21) 

which allows #13 to easily saturate its current upper limit. As another illustration, taking the current 
"best fit" values (6 12 , 23 ) = (34.5°, 42.8°) as in Table-1, we derive from (imp or (IPOip the lower 
limits 0i3 6° for the present model, and #13 ^ 1.5° for Ref. pQ. Hence, in contrast with Ref. p], the 
present model favors a larger $13, and can saturate its current upper limit, as will be demonstrated 
in Fig. 2 below. 

In the following, we systematically analyze the predicted parameter space and correlations in the 
present model (with inverted mass-ordering). We will find these to be very different from that in our 
soft breaking model (with normal mass-ordering) [1]. So, the present model can be tested against 
that in Ref. [1] by the on-going and upcoming neutrino experiments. 

Using the neutrino data for 6 S and (Am^, Am^ 3 ) (Table-1), and scanning the Dirac CP phase- 
angle 5 D G [0°, 360°) , we can plot the two fi— r breaking mixing angles, $13 (= 5 X ) and 023 — 45° (= 
5 a ), from (|4.12ap -( |4.12bp and (|4.13p . as functions of the theory parameter (^cosu and 8d • Our 
findings are depicted in Fig.QJa)-(d) with the experimental inputs varied within 90% C.L. ranges and 
with £cos uj G [—0.6, 0.6] in the natural perturbative region. Here we find that the theory prediction 
of 023 — 45° lies in the range, 

-4° < 023 - 45° ^ 4°, (4.22) 

which is within the current experimental bounds. On the other hand the predicted 0i3 can saturate 
the current experimental limits, and has distinct distributions. 
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Figure 1 : Predictions of #13 and 623 — 45° as functions of the fi — r breaking parameter ( cos uj 
and CP breaking parameter 5 D . The experimental inputs are scanned within 90%C.L. ranges and 
the Dirac phase angle 5 D E [0, 2ir) , with 1500 samples. The shaded region (yellow) denotes the 
90%C.L. limits on 8 13 and 23 - 45°, from Table-1. 

From the theory relations (|4.12a|) -( j4.12bp . we can further explore the correlation between the two 
fx — t breaking mixing angles #13 and 623 — 45° . This is displayed in Fig.[2j where we have varied the 
measured parameters within their 90% C.L. ranges, and input the Dirac-phase angle 5 D G [0, 2ir) as 
well as \('\ ^ 0.6 . The current 90% C.L. limits on 613 are shown by the shaded region (yellow), while 
the #13 sensitivities of the on-going Double Chooz [12] , RENO |T3] and Daya Bay [10] experiments are 
depicted by the three horizontal (red) lines at 90% C.L., as 5.0°, 4.1° and 2.9° (from top to bottom), 
based on three years of data-taking. The horizontal dashed (red) line represents Daya Bay's future 
sensitivity (2.15°) with six years of running [32]. 

Inspecting Fig. [21 we find that the sharp edges on the two sides of the allowed parameter space 
are essentially determined by the lower bound given in (|4.14p . S x ^ \5 a \ tan 29 s , where the current 
data require, 2.2 ^ tan2# s ^ 3.1 at 90% C.L. (Table-1) and the lower limit tan2# s = 2.2 just 
corresponds to the slopes of the sharp edges which are nearly straight lines. Hence, for any measured 
nonzero value of 623 — 45° 7^ , the Fig.\^ imposes a lower bound on O13 , which will be tested by the 
reactor experiments such as Daya Bay, RENO and Double Chooz. The current oscillation data favor 
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Figure 2: Correlation between #13 and #23 — 45° , based on Eqs. (|4.12ap -( |4.12b|) . where the experi- 
mental inputs are scanned within 90%C.L. ranges and the Dirac phase angle 5 D G [0°, 360°) , with 
1500 samples. The sensitivities of Double Chooz |12j . RENO [13] and Daya Bay |10j experiments to 
#13 are shown by the three horizontal (red) solid lines at 90% C.L., as 5.0°, 4.1° and 2.9° (from top 
to bottom). The Daya Bay's future sensitivity (2.15°) is shown by the horizontal dashed (red) line. 

the central value of #23 to be smaller than 45° (Table-1) and this feature is quite robust [6]. From 
Fig. [21 we see that taking the current central value of 823 — 45° = —2.2° (Table-1), the lower bound 
on #13 is already very close to the sensitivity of Double Chooz experiment; and a minor deviation 
of 023 — 45° = —1.4° will push #13 up to the sensitivity of Daya Bay experiment. Hence, the Daya 
Bay, RENO and Double Chooz reactor experiments hold great potential to discover a nonzero #13 . 
Furthermore, as shown in Fig.[2j detecting a nonzero #13 > 3° will strongly favor a nonzero #23 — 45° . 
Hence, we further encourage the improved measurements of #23 by Minos [7| and T2K [8], as well as 
future neutrino factory and super-beam facility |33} 134"] . 

Note that our previous soft breaking model [T] predicted a lower bound 8 X ^ \5 a \ tan# s with the 
slope 0.64 ^ tan# s ^ 0.73 at 90%C.L., which is about 3.4 — 4.2 times smaller than the present 
model. This means that given the same nonzero deviation of #23 — 45° , the current model will place 
a much stronger lower bound on #13 , higher than that in Ref. [I] by a factor of 3.4 — 4.2. Hence, the 
prediction of Fig. [2] is really encouraging for the upcoming neutrino oscillation experiments, which 
will probe the \i—t violating observables #13 — 0° and #23 — 45° to much higher precision. 

Then, we analyze our model predictions for the low energy CP-violation (via Jarlskog invariant J) 
and the neutrinoless double-beta decays (via the element |m ee | of M v ). From our theory construction 
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in Sec. 2. 2, the original CP-phase e lu) in the Dirac mass-matrix of seesaw Lagrangian is the common 
source of both low energy Dirac and Majorana CP-violations via the phase angles 5 D and 6<j> . 

The Dirac CP-violation is characterized by the Jarlskog invariant J [35] in the light neutrino 
sector with nonzero CP-phase 5 D and can be measured by the long baseline neutrino oscillation 
experiments. On the other hand, the neutrinoless double-beta decay observable \m ee \ contains both 
S D and Majorana CP-phase 5(f) . We can express the Jarlskog invariant J as follows [35], 

1 6 

J = - sin 20 s sin 26^8^26^ cos 6» x sin <5 D = -p sin20 s sm5 D + <D(8 2 X , 5 2 a ) , (4.23) 
8 4 

where as defined earlier, 5 X = O x and S a = 6 a — ? . The solutions (|4.12ap -( [4.12b[) leads to the 
correlation (|4. 13p . We can input the neutrino data for mixing angles, (9 S , 5 X ) = (#12, #13) , and 
mass-ratio, y' = Am^/ Am 2 3 , as well as scanning the model-parameter £' in its perturbative range 
IC'I < 0.6. 

We then study the neutrinoless double-beta decays. Our present model predicts the inverted 
mass-ordering (IMO) with m 3 = , so from (j3.5aj) we can derive the mass-matrix element \m ee \ for 
neutrinoless double-beta decays, 

M ee = \m ee \ =\52V e * j 2 m j = m l( ? x c 2 + s 2 y/l + y' e"* 2 * 



~ mi 



1 + \slv' -Si- 2s 2 s c 2 s l 



(4.24) 



where in the last step we have expanded 5 X and 5(p to the second order since y' = O(10 -2 ) is 
relatively small as constrained by the current data [cf. (|3.7p ], Eq. ()4.24p shows that the neutrinoless 
/3/3-decay observable M ee only contains the second orders of the fi—r breaking quantity 5 X (= 6*13) 
and the Majorana CP-phase angle 5cp . Hence, M ee is less sensitive to the fi — r breaking and 
Majorana CP-violation at low energies. 

We plot the correlation between 6*13 and the Jarlskog invariant J in Fig.[3ja), and the neutri- 
noless /3/3-decay observable M ee is depicted in Fig.[3jb). For the analysis of Fig.[3la), we have used 
Eq. (|4.12ap where we vary the model-parameter £' £ [—0.6, 0.6] in its perturbative range. We scan 
all other measured parameters within their 90%C.L. ranges. The shaded region (yellow) in Fig. [3] is 
allowed by the neutrino data at 90%C.L. Fig. E^a) shows that any nonzero J will lead to a lower 
bound on #13 due to 5 X 4| J|/sin26' s as inferred from Eq. (|4.23p . Combining the current upper 
limit 6*i3 < 9.5° (shaded region in yellow) with our parameter space in Fig.[3^b), we predict the 
allowed range, 

-0.037 < J < 0.037, (4.25a) 
45.5 meV < M ee < 50.8 meV, (4.25b) 

which can be probed by the on-going neutrinoless double beta decay experiments [3] . 

Before concluding this subsection, we compare our prediction (|4.13p with a recent independent 
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Figure 3: Correlations of #13 (in degree) with the Jarlskog invariant J [plot-(a)] and with the neu- 
trinoless /3/3-decay observable M ee [plot-(b)]. Each plot has computed 1500 samples. The shaded 
region (yellow) is allowed by the current data at 90% C.L. 



work [36]. In Ref. [36], using a charged lepton perturbation, Friedberg and Lee derived a very inter- 
esting prediction, cos 2#23 = tan 2 #13 , leading to 

1 



7T 
1 



7 23 



?? 3 « & 



13 



which does not contain CP phase and predicts a nearly maximal #23 
our predictions ([4.18ap -( |4.18b|) in the same notations, 



(4.26) 

For comparison, we rewrite 



7T 
1 
7T 
1 



703 



'23 



2 cos 5 d 



7T 



'12 



~ cos Sd 9 



13 



0(9 



13 



^13 < 6*13 , (Current Prediction), (4.27a) 
(Prediction of Ref. PQ), (4.27b) 



where our correlations explicitly contain the CP-phase angle 5 D . Moreover, our present model 
predicts a deviation j — 623 to be significantly smaller than #13 as in (|4.27aj) . due to the suppression 
of I — 9 12 = 0.16 — 0.21 at 90% C.L. But, taking cos<5 D = 0(1) , we see that the right-hand-side of 
(EOZaD is larger than that of by a factor of 4(f - 9 12 )/9 13 = (36.0 - 48.8°)/0i 3 at 90%C.L., 

which is clearly bigger than one. On the other hand, our previous soft breaking model [1] predicts 
the two small jjl— t breaking observabes to be of the same order, ^ — #23 = 0(^13) > as in (|4.27b|) . 
Hence, the predictions by Friedberg-Lee [36] and by us differ in a nontrivial and interesting way, 
which strongly motivate the on-going and future neutrino experiments for tests and resolution. 



4.2. Baryon Asymmetry from fx — r Blind Seesaw and Direct Link to Low Energy 

In this subsection, we study the predictions of our fi— r blind seesaw model for cosmological baryon 
asymmetry (matter-antimatter asymmetry) via thermal leptogenesis [37} [38] . We build up the direct 
link between leptogenesis CP-asymmetry and the low energy Dirac CP-phase, and further predict 
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the low energy leptonic Jarlskog invariant J [35]. Imposing the WAMP data on the baryon asym- 
metry [IB], we predict a negative Jarlskog invariant, J < , and derive a lower bound on the reactor 
mixing angle, #13 > 1° . We also analyze the correlations of the leptogenesis scale with the low en- 
ergy observables such as the Jarlskog-invariant J and neutrinoless /3/3-decay parameter M ee [3]. We 
further deduce a lower bound on the leptogenesis scale for producing the observed baryon asymmetry. 

Our universe is exclusively dominated by matter rather than antimatter. The asymmetry of 
baryon-anti-baryon density n B — n B (~ n B ) relative to the photon density n 7 is measured to be a 
tiny nonzero ratio [16] , 



r] B = — 2- = (6.19 ± 0.15) x 10" 10 . (4.28) 

The SM fails to generate the observed baryon asymmetry because of the too small CP-violations from 
CKM matrix and the lack of sufficiently strong first-order electroweak phase transition |39j , which 
violate Sakharov's condition for baryongenesis [3D] . It is important that the seesaw extension of 
the SM allows the thermal leptogenesis |37j with CP-violations originating from the neutrino sector 
and the lepton- number asymmetry produced during out-of-equilibrium decays of heavy Major ana 
neutrino Nj into the lepton-Higgs pair £H and its CP-conjugate IH* . Then, the lepton asymmetry 
can be partially converted to a baryon asymmetry via the nonperturbative electroweak sphaleron [31] 
interactions which violate B + L [32] but preserve B — L }43|, I44j . r] B = j-^b—l = ~\^L > wnere 
£ is the fraction of B — L asymmetry converted to baryon asymmetry via sphaleron process [43 1 and 
f = 28/79 for the SM. The dilution factor / = A 7 CC /A* = 2387/86 is computed by considering 
standard photon production from the onset of leptogenesis till recombination [33] . The effect of 
the heavier right-handed neutrino (A2) decays will be washed out in the thermal equilibrium, only 
the lightest one (Ai) can effectively generate the net lepton asymmetry for Mi <C M2 . (In the 
numerical analysis below, we will consider the parameter space with M%jM\ ^ 5 , to ensure the 
full washout of lepton asymmetry from A^-decays.) Thus, the net lepton asymmetry is deduced 
as [33], n[ = ^Kj-e 1 . Hence, we can derive the final baryon asymmetry, 



4/ 



Vb = — 7T K / e i = -dKf€ l , (4.29) 



where d = 3£/(4/) ~ 0.96xl0 -2 , and the factor Kj measures the efficiency of out-of-equilibrium Ai- 
decays. The is determined by solving the Boltzmann equation numerically [44|, 135] . In practice, 
useful analytical formulas for can be inferred by fitting the numerical solution of the Boltzmann 
equation. We find it convenient to use the following fitting formula of k^. [35], 

, / m, \ L16 3.3xl0~ 3 eV lA , 

with fn 1 = (rn} D rn£,) 11 /M± , and m D = m^U r with Ur being the rotation matrix diagonalizing 
the mass-matrix Mjj of right-handed neutrinos. In the present (jl — t blind seesaw, it is natural to 
set the right-handed neutrinos in their mass-eigenbasis from the start, Mr = diag(Mi, M2) , as we 
defined in Sec. 2. 1. So we have Ur = X with X the unit matrix, and thus m D = m D . (Other fitting 
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Figure 4: Leptogenesis scale M\ is plotted as a function of Dirac CP-phase angle 5 D , where the 
seven years of WMAP measurement (|4.28p is imposed. All experimental inputs are scanned within 
their 90% C.L. ranges, with 1500 samples. 



formulas than (|4.30j) to the exact solution of Kj in the literature [H] agree with each other quite 
well for the relevant range of m 1 .) The CP asymmetry parameter e 1 is defined as 

rjJVi -T\Ni -+IH*\ 1 F (M 2 \ 9m {fefl)i2l 2 } 

where v denotes the vacuum expectation value of the SM Higgs boson. As we constructed in Sec. 2. 2, 
the Dirac mass-matrix m D is complex and provides the common origin of the //— r and CP breaking; 
the complexity of m D causes the difference between the decay widths T[N\ —> £H] and T[N\ — > 
£H*\ , and thus a nonzero CP asymmetry e x 7^ . For the SM, the function F(x) in (|4.3ip takes the 
form, 



F{x) = x 



1 + x 2 1 



l-(l + x 2 )ln ^ + 



x 2 1 — x 2 



For our numerical analysis of the thermal leptogenesis, the mass ratio M2/-M1 3> 1 and thus the 
above expanded formula of F(x) holds with good accuracy. 
Then, we proceed to compute the matrix elements, 

(m\ ) m D ) 11 = m Mi (a 2 + 2b 2 ) = m Mi , (4.33a) 

{m j D m D ) 12 = -mo^/MjhbciC + (e iuj ) . (4.33b) 

So we can deduce the effective mass-parameter m 1 as introduced below ([4.301) . 



l 13 ' 

and the imaginary part, 



mi = fh ~ Xi\/Amf 3 , (4.34) 



Sm|[(m^m D ) 12 ] 2 } = ~ rh 2 M x M 2 y' sm29 s sin<5 D S x , (4.35) 
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Figure 5: Seesaw scale M\ and M 2 as functions of the elements (a, 5) and (a', c) in the Dirac mass- 
matrix m D , where the shaded regions correspond to the natural perturbative region (a, b, a', c) £ 
[1, 300] GeV, and 600 samples are generated in each plot. This puts an upper bound, M\ ^ 3.5 x 
10 15 GeV from plot-(b), and M 2 < 1.7 x 10 15 GeV from plot-(c). 



where the RG running factor \\ = x(Mi,m z ) is defined in Eqs. (|3.18p - (|3.19p . Using Eq. (|4.34p 
together with the neutrino data (Table-1), we find that the light neutrino mass-parameter m 1 lies 
in the 3a range, 0.046 < m l /x\ < 0.053 eV, where the RG factor Xi — 1-3 ~~ 1-4 is evaluated 
numerically, as explained around the end of Sec. 3. 2. So, in Eq. (|4.30p the second term on the 
right-hand-side is negligible and Kj is thus dominated by the first term. 

With these and from (|4,3ip . we derive the CP asymmetry parameter e 1 as follows, 

3y'rh Mi 



ci - 16th; 2 sin26>s sinS D $r 

Finally, inspecting Eqs. ()4.29p . (|4.3ip and (|4.32p . we can derive, 



sin 29 s sin 5 D 5 X 



(4.36) 



(4.37) 



Mi ' vf 167TV 2 

Since the WMAP measurement (|4.28p finds the baryon asymmetry ij B > , so we can infer the 
constraint, sin 5 D < 0, which restricts the Dirac phase angle, 5 D £ (ir, 2ir) . 

Then, from Eq. (|4.37p we compute the ratio T] B jM\ for any nonzero sin^ , where we vary all 
measured quantities within their 90% C.L. ranges. Since < | sin<5 D | ^ 1 , we can deduce a robust 
numerical upper bound, 



-?f < 1.8 x HT^GeV- 1 . 
Mi 



(4.38) 
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Inspecting (|4,37p we can also reexpress the leptogenesis scale M\ in terms of baryon asymmetry r] B 
and other physical observables, 

jyj^ _ 39^ 

3d Kjrin y' sin 20 s sin Sd S x 

With the data of r] B from (|4.28[> . we can plot, in Fig.[U the leptogenesis scale Mi as a function of 
Dirac CP-phase 5q , where all experimentally measured quantities are scanned within their 90% C.L. 
range (with 1500 samples). Fig. H] reveals a robust lower bound on Mi , 

Mi > 3.5 x 10 13 GeV . (4.40) 

Using Eqs. ()2.7|) and (|4.6p . we connect the seesaw scale (Mi, M2) to the elements of the Dirac 
mass-matrix m D , 

a 2 2 b 2 

Mi = — tt— — = = , (4.41a) 

mo cos z 20 s mo sin 26 s 

a' 2 2 c 2 

M 2 = - % = - Z \ , (4.41b) 
mo sin 20 s uiq cos z 20 s 

where the Dirac mass-parameters (a, b, a', c) arise from the Yukawa interactions, (a, b, a' , c) = 

(?/a; 2/a'' yc) v /^- So we can plot Mi as a function of the magnitude of the Dirac mass-parameter 

\a\ or |6| in Fig.[S]^a)-(b), and M2 as a function of the magnitude of the Dirac mass-parameter \a'\ 

or |c| in Fig.[5jc)-(d), where we have varied the measured quantities in their 90% C.L. ranges. We 

note that the Yukawa couplings (y a , y b , y a ,, y c ) cannot be too small (to avoid excessive fine-tuning) 

or too large (to keep valid perturbation). So, we will take the Dirac mass-parameters (a, b, a' , c) in 

the natural range [1, 300] GeV, corresponding to the Yukawa couplings y^ no smaller than O(10~ 2 ) 

and no larger than 0(y t ), where y t = \f2m t jv ~ 1 is the top-quark Yukawa coupling in the SM. 

This natural perturbative range of (a, b, a', c) is indicated by the shaded area in Fig.[5ja)-(d), which 

results in an upper limit on the seesaw scale (Mi, M2) due to the perturbativity requirement. From 

Fig.[5^b) we infer an upper bound Mi < 3.5 x 10 15 GeV, while Fig.[5^c) requires M2 1.67 x 

10 15 GeV. For the above construction of natural thermal Leptogenesis we consider the parameters 

space M2/M1 5 , so with the upper bound of Fig.JS^c) we further deduce a stronger limit Mi ^ 

3.3 x 10 14 GeV. 

With the above constraint on the parameter space from realizing successful thermal leptogenesis, 
we can rederive the correlation between #13 and #23 — 45° , as shown in the new Fig.[6l which should 
be compared with Fig. [2] in Sec. 4. 1 (without requiring leptogensis). We note that the realization of 
successful thermal leptogenesis puts a general lower bound on the mixing angle #13 , 



'13 



> 1° , (4.42) 



even for the region around #23 = 45° . 

Under successful leptogenesis, the correlations of #13 with the Jarlskog invariant J and the 
neutrinoless double beta decay observable M ee are plotted in Fig.[7Ja) and (b), respectively. This 
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Figure 6: Correlation between #13 and #23 — 45°, where all the inputs are the same as Fig.[2l except 
requiring successful leptogenesis in the present analysis, with 1500 samples. 

should be compared to Fig. [3] where leptogenesis is not required. We see that due to the constraint 
from the observed baryon asymmetry, the parameter space of J > is forbidden in Fig.[7^a). On the 
other hand, the constrained range for M ee in Fig.[7^b) is almost the same as Fig.[3^b), since Eq. (|4.24|) 
shows that the observable M ee has rather weak dependence on small NLO parameters S x (= #13) 
and 5(j) via their squared terms. Thus, from Fig.[3a)-(b), we infer the following constraints on J 
and M ee , 



which should be compared to Eqs. (|4.25ap -( |4.25b|) in Sec. 4. 1 without requiring the successful lepto- 
gensis. 

We further analyze the correlations of the neutrinoless /3/3-decay observable M ee with the Jarlskog 
invariant J and the light neutrino mass m 1 (~ m 2 ) , in Fig.[8^a-b) and Fig.[8](c-d), respectively. The 
two left plots in Fig.[8ja) and (c) show the correlations of M ee with J and with m 1 after imposing the 
leptogenesis. For the two right plots in Fig.[5^b)(d), we have replotted the same model-predictions 
as in the two corresponding left plots of Fig.[8^a)(c) (all in blue color). For comparison, we have 
further plotted, in Fig.[8^b)(d) with green color, the model-independent parameter space of M ee 
[cf. (|4.24p ] versus J [cf. (j4.23|) ] or m 1 (= \J Amf 3 ) , for the IMO scheme with 7713 ~ , where the 
relevant observables are varied within their 90%C.L. ranges and 5 D G (0, 2ir]. This comparison 



0.037 < J < - 0.0035 



(4.43a) 



45.5 meV < M ee < 50.7meV 



(4.43b) 
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Figure 7: Correlations of #13 with Jarlskog invariant J in plot- (a) and with neutrinoless double 
beta decay observable M ee in plot-(b), where all inputs are the same as Fig.O except requiring the 
successful leptogenesis in the present figure, with 1500 samples for each plot. 



shows that our model predictions are located at the upper boundaries of the whole parameter space, 
giving rise to the largest allowed M ee . This is very distinctive and highly testable. Furthermore, in 
Fig.[5](b)(d), we have compared our predictions with the sensitivities of the future neutrinoless j3j3- 
decay experiments CUORE (CU) @Bj and Majorana @7j/GERDA III 08] (M/G), which are depicted 
by the horizontal dashed lines at 15meV (black) and 20meV (red), respectively. 

The leptogenesis scale M\ can be determined from the baryon asymmetry r] B , the reactor angle 
#13 , the Dirac phase sin 5 D and other neutrino observables as in Eq. (|4.39p . Since the low energy 
parameter J in Eq. (|4.23p is also predicted as a function of #13 and sin 5 D , so it will correlate 
with the leptogenesis scale M\. Hence, we can plot the correlations of the leptogenesis scale M\ 
with the reactor angle #13 in Fig.[9^a), and with the Jarlskog invariant J in Fig.[9^b). Inspecting 
Eqs. (|4.23p and (|4.39|) . we deduce, J cx b~ x smbD and M\ cx (5 X sin^o) -1 , from which we arrive 
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Figure 8: Upper plots (a)-(b) show the correlations between the neutrinoless /3/3-decay observable M ee 
and the Jarlskog invariant J with successful leptogenesis. Lower plots (c)-(d) depict the correlations 
between M ee and light neutrino mass m 1 (~ m 2 ) with successful leptogenesis. All experimental 
inputs are varied within 90% C.L. ranges, for 1500 samples. The background (green) regions in plots 
(b) and (d) represent the model-independent parameter space of the IMO scheme with 7713 ~ . 
The horizontal dashed lines in (b) and (d) depict the sensitivities of the future neutrinoless /3/3-decay 
experiments CUORE (CU) [IB] and Majorana 0Z1/GERDA III 08] (M/G), at 15meV (black) and 
20meV (red), respectively. 



at, Mi oc 1/1 J I . This behavior is impressively reflected in Fig.[9jb), as expected. In addition, the 
relation, M\ oc (S x sin<5z)) -1 ^ #13 , nicely explains the lower arched edge in Fig.[9^a). 

4.3. Extension to General Three-Neutrino Seesaw 

In this subsection, we analyze the extension to the general neutrino seesaw with three right-handed 
neutrinos j\f' = (N\, N2, N^) T , where M' is fi— r blind. Then, in the ll—t and CP symmetric limit, 
the mass-matrices m D and Mr are extended to 3 x 3 matrices, 



m 



D 



la a' 


a") 




( a i a 


a 2 d 


a 3 a" 
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b c 


d 




a^b 


a 2 c 


a 3 d 




\b c 


d j 




\ *ib 


a 2 c 


a 3 d 


) 


^fh M 1 


, a 2 


= y/m Q M 2 


, and 


^3 = 





M R = diag(Mi, M 2 , M 3 ), (4.44) 



neutrinos N 1 can always be rotated into their mass-eigenbasis without affecting the structure of 
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Figure 9: Correlations of leptogenesis scale Mi with the reactor mixing angle #13 in plot-(a), and 
with the low energy Jarlskog invariant J in plot-(b). Each plot contains 1500 samples. 



m D . Thus, we rederive the \i — r and CP symmetric seesaw mass-matrix for the light neutrinos, 
(a 2 +a' 2 +a" 2 ab + a'c + a"d ab + a'c + a"d\ / A B B 



Win 



b 2 +c 2 +d 2 



V 



b 2 +c 2 +d 2 
b 2 +c 2 +d 2 J 

from which we deduce the mass-eigenvalues and mixing angles, 

1 



C s C s 



(4.45) 



m 



1,2 



mo 
2 



(A + 2C S ) =f ^{A - 2C S ) 2 + 8B 2 
[(a 2 + a' 2 + a" 2 + 2b 2 + 2c 2 + 2d 2 



=F V[(« 2 + a' 2 + a" 2 ) - 2{b 2 + c 2 + d 2 )] 2 + 8{ab + a'c + a"df 



tan 2^i2 



Cs ~ C s — , 
2y/2B s 



2 v / 2|a& + a'c + a"d\ 



A-2C S ' \a 2 +a' 2 +a" 2 -2{b 2 +c 2 +d 2 )\' 
45° , 013 = 0° , 



(4.46a) 
(4.46b) 
(4.46c) 
(4.46d) 



where the mass-spectrum remains the inverted mass-ordering (IMO). The third mass-eigenvalue m 3 
vanishes because our \i — r blind seesaw (|4.44|) predicts the seesaw mass-matrix (|4.45|) with its 23- 
element equal to the 22-element and 33-element. This is also a general feature of any fi—r symmetric 
IMO scheme at the LO, as to be shown in (|5.5p of Sec. 5.1. Furthermore, we will demonstrate shortly 
that the third mass-eigenvalue m 3 = actually holds up to the NLO after including the fi— r and 
CP breaking in our analysis. So this resembles very much the minimal seesaw we studied earlier. 

Similar to Eqs. (|2.1ip and (|2.13|) in Sec. 2. 1, we can realize the IMO at the LO of three-neutrino 
seesaw, m 1 = m 2 = m , which leads to the three extended conditions, 



( a 2 + a '2 + fl "2) + 2 (6 2 + c 2 + d 2 } = 2) 
( a 2 + a '2 + a //2 ) _ 2(6 2 + c 2 + d 2 ) = Q) 



(4.47a) 
(4.47b) 
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ab + a'c + a"d = 0. 



(4.47c) 



With these we deduce from (|4.45p the generic LO seesaw mass- matrix for the IMO, 

/ a 2 +a' 2 + a" 2 ab+a'c+a"d ab+a'c+a' 



lllr 



b 2 + c 2 + d 2 



V 



b 2 + c 2 +d 2 

u2 , „2 , j2 
















i 








2 


i) 




I 




2 / 



, (4-48) 



which is the same as the LO mass-matrix (|2.31ap we derived earlier for the minimal seesaw. Hence, 
despite that the LO mass-matrix Mo contains two new parameters (a", d) at the beginning, the 
realization of IMO eliminates them all and reduces M„ to the universal LO mass-matrix as shown 
in the final form of (|4.48p which is parameter-free except an overall mass-scale. As a result of the 
IMO conditions (|4.47|> . we note that the solar angle formula (|4.46cp gives tan2#i2 = j at the LO, 
which is now undetermined. So, the 812 has to be derived from the NLO contributions related to 
li — t breaking terms. Before getting into detail, it is convenient to infer 9\2 by using the l'Hopital 
rule, similar to what we did in Sec. 2. 1 for the minimal seesaw. Thus we have, 



tan 29 



12 



y/2\b\ ' 

for li — t breaking arising from the deviation in the element b of m D , or 

\n'\ 

tan 2#i2 



(4.49a) 



V2\c\ ' 

for li — t breaking arising from the deviation in the element c of m D , or 



(4.49b) 



tan 29 



12 



V2\d\ 



(4.49c) 



for li — t breaking arising from the deviation in the element d of m D . 

As noted in Sec. 2.2, we can always rotate the first column in m D to be all real by rephasing. For 
the convenience of comparison with the minimal neutrino seesaw, we will thus formulate the common 
origin of li—t and CP breaking in the element c of m D . It is possible to construct such a breaking 
in the element d of m D , but this does not affect our physical conclusions as will be clarified below, 
after Eq. (|4.57p . [Since we are constructing a common origin of li — t and CP breaking from a single 
source in m D , we do not consider this breaking to occur in both c and d elements of m D at the 
same time.] So, we build the Dirac mass-matrix m D with the common li—t and CP breaking in the 
following form, 

/ 

m D = 





a 2 a' 


a,a" \ 




°2 C 1 


a 3 d 




02 c 2 


a 3 d J 


-0 


c 2 


= c(l- 



Ce-) 



(4.50a) 
(4.50b) 
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Thus we can deduce the NLO part of the seesaw mass-matrix M v = M„ -\-5M„ for light neutrinos, 



6MP 



mr 



/0 -a'cC -a'c(e iuj \ 

-2c 2 C -c 2 {(' + (e iui 
y -2c 2 (e iw J 



(4.51) 



which equals (|2.31b|) as expected, since the new parameters (a" , d) appear in the seesaw mass-matrix 
M u only via the products (a" 2 , d 2 , a"d) with no crossing terms like c x 2 a" or c 1 2 d- With these, we 
deduce the \i—t symmetric and antisymmetric elements of 5M^ to be the same as Eq. (|2.33p . 

Using the formalism of Sec. 3. 1 and extending Sec. 3. 2, we can reconstruct the light neutrino 
mass-matrix M v for the IMO with m 3 7^ , via the NLO parameters, 



(y, z, z, 5 a , 5 X , 5ai, 5a 2 , 5a 3 , 



0, 



(4.52) 



where we have defined z' 

1 / / 



in.j 



in. 



and 



!>i = </>n + S<f>' . Note that the LO phases vanish, 



a 



3 — ^1 — ^0 

— . So the NLO elements of M u are reconstructed as follows, 

2 

. / „-r»/_2 ri , Cr 

I, 



5A 



8B S 

SC S 

5D 
SB a 

5C a 



z + yy' - i2(s 2 5<j> + 6 



2V2 



Wr 



sin 29 s 



1 , . 

-y + % 
2 y 



z + z' + ^y'- i(2c 2 



z-z - i{2c s 



+ 5a 2 + Sets) 
+ 8a 2 + 5a 3 ) 



V2 X 
-mo 



5a + ^(^a - 8CX3) 



(4.53a) 

(4.53b) 

(4.53c) 

(4.53d) 
(4.53e) 

(4.53f) 



where we note that the Majorana phase 5(p' does not appear at the NLO because it is always 
suppressed by another NLO parameter z' = ^f- . Moreover, since the \i— r and CP breaking matrix 
(|4.5ip gives Eq. (|2.33|) with the equality 5C S = 5D , we deduce z' = — a = by comparing (|4,53c-p 



with ()4.53dp . and thus 



ITln 



holds up to the NLO. Hence, we have shown that our model with 



the general three-neutrino seesaw under IMO does share the essential feature of m 3 = with the 
minimal seesaw. 

Then, with the NLO fx— r symmetric parts from (|2.33p and (|4.53p . we deduce the solar angle #12, 



tan 29* 



V2c ' 

which coincides with Eq. (|4.5p as we derived earlier for the minimal seesaw. 
Next, connecting the fj,— r anti-symmetric parts in (|2.33|) and (|4.53|) gives, 



(4.54) 



L r iSp 

V2 



e~"6 x , 



(4.55a) 
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from which we arrive at 



m c 2 (C'-Cen 

cos 5d 8 x 

sin 5 D 5 X 
S a 

6ct2 — 8a% 



-mn 



a'c 



S a + 2(^*2 - 80*3) 



-^■(C'-Ccosw) 



a'c 



7f ^ sina; ) ' 



■ 2 (C'-Ccosw), 
-2c 2 (C sin a;) . 



(4.55b) 

(4.56a) 

(4.56b) 

(4.56c) 
(4.56d) 



Here for the left-hand-sides of ()4.55ap -( |4.55bj) we have used the Eq. ()3.19p to evolve the overall mass- 
parameter m from seesaw scale down to the corresponding m at low energy. 

Finally, using Eqs. (|4.54|) . (|4.56a|) and (|4.56c|) . we derive the key correlation between two low 
energy (j,— r breaking observables S a and 8 X , 



5 a 



cot 26 s cos 5d 8 x , 



(4.57) 



which coincides with ()4.13p as we derived earlier for the minimal seesaw. 

We note that it is also possible to construct the common origin of fi— r and CP breaking in the 
element d of m D , instead of the element c . Then we can rewrite the Dirac mass-matrix ()4.50|) as 





( 




o 2 d 


a,a" 


) 




m D = 






a 2 c 


Ml 




(4.58a) 




{ 




a 2 c 


a 3 d 2 


J 





d 1 = d{\- CO , d 2 = d{\- Ce luJ ) 



(4.58b) 



This results in the following NLO seesaw mass-matrix, 

/0 -a"d(' 

8M^ = m 



\ 



a"d(e iw \ 
2d 2 C -d 2 (('+(e iui 



(4.59) 



-2d z (e tLU J 



from which we derive the solar angle, 



tan 20, 



V2d ' 



and the reconstruction conditions, 



COS Of) 5 X 

sin5 D 5 X 

8a 



a"d (H \ 
(C -Ccosw) 



a"d 



(Csinw) , 



V2 

d 2 (C'-Ccosw) 



(4.60) 

(4.61a) 

(4.61b) 
(4.61c) 
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<5o2 — Sas = — 2d 2 ((, sin uj) 



(4.61d) 



So, from Eqs. (|4.60p . ()4.61ap and ()4.61cp . we can readily derive the correlation between two [i — r 
breaking observables, 

5 a = — cot 26 s cos 5d S x , (4.62) 

which coincides with (|4.57p . 

In summary, the general three-neutrino seesaw (with right-handed neutrinos being [i — t blind) 
still predicts the inverted mass-ordering (IMO) for light neutrinos [cf. Eqs. (|4.46ap -( |4.46bp ]. Despite 
that the LO conditions (|4.47|) for the IMO contains two new parameters (a", d), the LO seesaw 
mass-matrix (|4.48p is shown to take the same form as in the minimal seesaw. Furthermore, the NLO 
[i — t and CP breaking part of our seesaw mass-matrix (|4.5ip or ()4.59p exhibits the same structure 
as in the minimal seesaw. This makes our final physical prediction of the correlation (|4.57p or (|4,62p 
coincides with (|4.13p . 

5. Hidden Symmetry and Dictation of Solar Mixing Angle 

So far, by analyzing the fi— r symmetry and its breaking, we have studied the atmospheric mixing 
angle 623 and the reactor mixing angle #13 in great detail. As shown in Table-1, the solar mixing 
angle 612 is best measured |49|. [50] among the three mixing angles. In this section we will clarify 
the connection between fj,— r breaking and the determination of the solar mixing angle 612 for both 
inverted mass-ordering (IMO) (cf. Sec. 2) and normal mass-ordering (NMO) [T]. Then, we analyze 
the general model-independent Z2 ® 1>2 symmetry structure in the light neutrino sector, and map it 
into the seesaw sector, where one of the Z2 symmetries corresponds to the fi— r symmetry and 
another the hidden symmetry Z| (which we revealed in [1] for the NMO of light neutrinos and is 
supposed to dictate 612 ). We will further derive the general consequences of this Z| and its possible 
violation in the presence of \i—t breaking for cases either with or without neutrino seesaw, regarding 
the d\2 determination. 

5.1. \i — t Breaking versus #12 Determination: Inverted Mass-Ordering 

In Ref. [1] we proved that the solar mixing angle #12 (= 9 S ) is not affected by the soft [i—t breaking 
from the neutrino seesaw, and we revealed a hidden symmetry Z| for both the seesaw Lagrangian 
and the light neutrino mass-matrix which dictates 9 S , where the normal mass-ordering (NMO) is 
realized. In this subsection, we generally analyze mass-eigenvalues and mixing angles for the jjL— r 
symmetric mass-matrix of light neutrinos under the inverted mass-ordering (IMO). Then we explain 
why the fi—r breaking is invoked for the 9 S determination and why the hidden symmetry Z| will be 
violated. The /i — r blind seesaw constructed in Sec. 2 belongs to an explicit realization of the IMO 
scheme. 
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Let us start with the general fi—r symmetric mass-matrix for light neutrinos, 

A B x B x 



which can be diagonalized as follows [TJ [S3] , 



C s D 



(5.1a) 



1.2 



tan 2$, 



[A + (C s +D)} T\/[A- (C s +D)} 2 + 8S2 

C s -D, 
2y/2B t 



A-(C S +D) 
Substituting ()5.2cp into ()5.2ap . we arrive at 



45° , e a 



0° 



??? 



1,2 



i {[A + (C s +L>)] T I A - (C fl +D)| sec20 s } 



(5.2a) 
(5.2b) 

(5.2c) 
(5.3) 



For the IMO scheme, we have the mass-spectrum m 2 > m 1 S> m 3 , where a small m 3 ^ is also 
generally allowed for the analysis below. So we can derive, for the general IMO scheme, 

A- (C s +D) 



A + (C s +D) 



m — mi 



cos 26 s 



Aran, 



(2.1 - 3.8) x 10" 



(5.4) 



m 2 + m 1 

where in the last step we have used the neutrino data (Table-1) to estimate the allowed range of 
this ratio at 90%C.L. Literally, Eq. (|5.4|) shows a fine-tuned cancellation between the mass-matrix 
elements A and (C s + D) down to the level of 10~ 3 . As will be clear in Sec. 5.2. 2 by using the 
general reconstruction formalism for the IMO scheme, we find that the LO form of the fi— r symmetric 
mass-matrix predicts the exact relations [cf. Eq. (|5.44p ]. 



#)-(Cf) + #) = 
J3(°)=0, Cf)-L>(°) 



1 1 

2 + 2 



0. 



(5.5) 



0. 



which ensures m 1 = m 2 and ?n 3 = at the LO. So, the small ratio (|5.4p naturally arises from 
the NLO elements [5 A — (5C S + SD)] ^ , and thus there is no real fine-tuning in (|5.4p . This also 
means that at the LO the solar angle 9 S is undetermined from the formula (|5.2cp , tan 29 s = S , and 

(s) 

the real determination of 6 S is given by the NLO elements of Mi, , 

2\/2 8B S 



tan 26* 



(5.6) 



6A- (5C S +5D) ' 

as we will explicitly verify in the next subsection for the general IMO scheme [cf. Eqs. (|5.46|) - (|5.47a|) ]. 

For the fi— r blind seesaw defined in Sec. 2. 1, we find that the light neutrino mass-spectrum must 
be inverted ordering, as given in Eqs. (|2. 10a.|) -f |2.10b[) . So, following the consistency with neutrino 
data (|5.4p and matching the reconstruction formalism (|5.5p for the IMO scheme, we can explicitly 
realize the degeneracy m l = m 2 at the LO by imposing the condition ()2.1ip on the elements of m D . 
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(Here m 3 = is an outcome of the minimal seesaw.) Thus, as expected, we find a problem for the 
9 S determination in the fi — r symmetric limit, 

tan2#s = J j ^— " = 7; > (5-7) 

which is just an explicit realization of our above general IMO analysis [cf. (|5.5|) ]. Hence, it is clear 
that 9 S must be inferred from the NLO formula (|5.6p . where the NLO elements will be predicted 
by a given model, e.g., by the first four expressions in Eq. (|2.33p in the fx— r blind seesaw with all 
NLO corrections arising from the fx— r breaking [51] . Thus the explicit expression of 9 S from such 
underlying models will depend on how the fx — r breaking is constructed. This is contrary to the 
neutrino seesaw with normal mass-ordering (NMO) of light neutrinos as studied in Ref. pQ, where we 
find that the formula of tan2# s [cf. (|5.2cj) above] is well defined in the fx— r symmetric limit. 

As we noted in Sec. 2. 1, the structure jj in Eq. (|5.7p allows us to use the l'Hopital rule on (|5.7p 
by taking the first derivatives on both its numerator and denominator. We need to decide for which 
parameter in (|5.7p the derivatives should be taken. There are only two possible choices, either c or 
b, since the fx— r breaking under the fx— r blind seesaw could appear in either c or b element of m D , 
as we explicitly constructed in Eqs. f|2. 18|) and (|2.23p . Thus, applying the l'Hopital rule to (|5.7p we 
have 

■ , (fx— t breaking in c) , 
tan2# s = { , C ' (5.8) 



a 



(fx — t breaking in b) , 



V2\b\ ' 

which, as expected, gives finite expressions for 9 S , depending only on the LO parameters of the 
Dirac mass-matrix m D . This also agrees to Eqs. (|2.15p - (|2.16p in Sec. 2. 1. But Eq. (|5.8|) shows that 
9 S does depend on how the fx— r breaking is built in the seesaw Lagrangian, and the two different 
constructions of fx— r breaking for m D lead to two different 9 S formulas above. This is an essential 
difference from the soft li—t breaking model in Ref. [1], where 9 S is dictated by the hidden symmetry 
Z| under which the soft fi— r breaking term in Mr is an exact singlet. In the next subsections we 
will analyze the general model-independent Z2 <8> ^2 symmetry in the light neutrino sector, and then 
map it into the seesaw sector. This allows us to explore, at a deeper level, the Z| symmetry and its 
possible partial violation under the fx — T breaking in a unified way, concerning 9 S determination. 

5.2. Z2 Symmetry under General fx — r Breaking and General Determination of 

This subsection consists of two parts. In Sec. 5. 2. 1, we analyze the general model-independent 
symmetry structure of the light neutrino sector, in both the mass and flavor eigenbases. We will show 
that, in the flavor eigenbasis of light neutrinos, one of the Z2's is the U£ symmetry which predicts 
the mixing angles (#23; $13) = (45°, 0°) , and another is the Z| symmetry which generally dictates 
the solar angle 9 12 by its group parameter (allowing deviations from the conventional tri-bimaximal 
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mixing ansatz). With general fJL — T breaking parameters, we will derive a nontrivial correlation 
between the two fi — r breaking observables which is necessary for holding the Z| symmetry. In 
Sec. 5. 2. 2, we will further analyze the general fj, — r breaking in the light neutrino mass-matrix M v 
and derive a nontrivial consistency condition to hold the Z| symmetry. From this condition and 
using the general reconstruction formalism of Sec. 3. 1, we will deduce the same correlation between 
the \i— t breaking observables, for both the normal mass-ordering and inverted mass-ordering of light 
neutrinos {without approximating the lightest neutrino mass to zero) [52]. 

5.2.1. Z| Symmetry for General Determination of Solar Angle #12 

Let us inspect the flavor symmetries in the lepton and neutrino sectors. In general, the lepton 
and neutrino sectors are expected to obey different flavor symmetries. After spontaneous symmetry 
breaking, the residual symmetry groups for the lepton and neutrino mass-matrices may be denoted 
as Gi and G v , respectively. Consider the symmetry transformations Fj 6 G^ and Gj £ G v for 
left-handed leptons and neutrinos. Thus the mass-matrices of leptons (Mi) and light neutrinos (M u ) 
will satisfy the invariance equations [57] , 

FjM e MjFj = M(M[ , GjM u Gj = M v . (5.9) 

The above mass-matrices can be diagonalized by unitary rotations for left-handed leptons and neu- 
trinos, 

u\M t M\Ui! = D e = diag (mf, m^, m 2 T ) , UjM u U u = D u = diag(m 1 , m 2 , m 3 ) . (5.10) 

Then, combining the invariance equations (|5.9p and diagonalization equations (|5.10p result in 

U}F}M t M}F j U e = d\D e d e = D e , U^G] ^M u G j U u = a%D v d v = D u , (5.11) 

where d £ and d u are diagonal phase-matrices obeying d\d^ = Z3 and d\ = X 3 (with X 3 the 3x3 
unit matrix), which require dp = diag(e* 71 , e* 72 , e* 7 3) and d u = diag(±l, ±1, ±1) . So, up to an 
overall phase factor, the {gQ j ^} forms the generic Abelian group U(l) ®U(1) = Gp for leptons, and 
{dv } has only two independent d v , 

d« = diag(l, 1, -1) , 4 2 ) = diag(-l, 1, 1) , (5.12) 

forming the generic discrete group Z2 Z2 — Gp for neutrinos. From Eq. (j5. 1 ip the following 
consistency solutions are deduced, 

Fj = U e d { e j) uj, G 3 = U v d^Ut. (5.13) 

This proves that {Fj} and {d^} are just connected by the similarity transformations, and are thus 
two equivalent representations of the same group ; similarly, {Gj} and {d^} are two equivalent 
representations of the same group G u . We may call the representation {d^p} and {dP} the "kernel 
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representations", with which the equivalent "flavor representations" {Fj} and {Gj} can be generated 
as in (|5.13|) via the disgonalization matrices U% and U u , respectively. Hence, we are free to choose 
an equivalent lepton symmetry group representation {Fj} = {dg } with U e =1 3 , and accordingly, 
rewrite the representation of neutrinos symmetry group, 



G ; 



vdi j) W, 



(5.14) 



with V = U\ U u = U u equal to the physical PMNS mixing matrix as defined in Eq. f|3.3|) of Sec. 3[ 
Let us rewrite the PMNS matrix (|53|> . V = U"UU' = V'U' , with V = U"U as introduced in 
Eq. (|3.4p . So we see that the Majorana phase-matrix U' cancels in Gj, 



G, 



v'dp 



(5.15) 



According to the most general reconstruction formulation in Sec. 3. 1, we can expand the matrix 
V' to NLO in terms of the small parameters, (5 a , 5 X , Sctj) , where (S a , 5 X ) characterizes the low 
energy fj,— r breaking and the CP-angle 5a i arises from the phase matrix U" (which is not directly 
observable and only needed for the consistency of diagonalizing the mass matrix M u ). There is no 
need to expand the Dirac CP-phase e lS ° itself since it is always associated with the small // — r 
breaking parameter 5 X . So, under this expansion we derive 



V' = V S + 8V'. 



(5.16) 



with 



V* 



1 c s 


-s s 


o \ 


Ss 


Cs 
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V2 


V2 




Ss 
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Cs 


7f / 


V2 



(5.17a) 
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ic s 5oti 

SsS a +CsS x e lS D +is s 5ct2 

s s S a +c s S x e' s D -tSsSas 
V2 



is s 8a\ 

-CsSa+s s S x e tS D -ics8a.2 
CsS a -s s S x e t6 D -iCsSas 



V2 



-6 x e- iS »\ 

8g—i8a2 

_ Sg+i&as 
V2 / 



(5.17b) 



Let us first consider the fi— r symmetric limit with V = V s . So substituting V s into Eq. (|5.15p we 
deduce, 



Gjj, T = G 1 




(5.18) 



which, as expected, just gives the Zg T symmetry-transformation matrix G v for light neutrinos as we 
explicitly constructed in (|2.5p earlier for the seesaw Lagrangian (|2.2|) . 

Next, we derive the symmetry-transformation matrix G® corresponding to of (j5. 12j) in the 
[i — t symmetric limit with (5V = 0), 

-y/2s s c s -V2s s c s \ 



s 2 -c 2 



Gri 



(5.19a) 



/ 
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1 + k 2 



which is symmetric since V s and d„ are real. (For the same reason G^ T is also symmetric.) In the 
last step, for convenience we have defined, 



(k 2 -\ 



\ 



-V2k 
1 



-V2k\ 
-k 2 
1 / 



(5.19b) 



^ s s> VT+ti 2 

with k (or equivalently, tan# s ) serving as the group parameter of 



(5.20) 



J 2 > 



tan# s 



k. 



where we can always choose the convention of 9 S G [0, ^] such that, tan# s = 
(du^) 2 = I3 and using the relation G° = V s d^Vs , we can readily verify (G®) 2 
deed G® £ 7L% . Hence, the solar angle 6 S is dictated by the group parameter k of the 3-dimensional 



(5.21) 

k ^ 0. Noting 
= X3 and thus in- 



representation of the hidden symmetry Z| [55] . We stress that the G® in (j5. 19bj) , as the 3d repre- 
sentation of Z| , is uniquely fixed by the fj,— r symmetric matrix V s ; we call Z| a hidden symmetry 
since it generally exists for any /U — r symmetric neutrino mass-matrix Mu [cf. Eq. (|5.29ap below], 
i.e., any \x—t symmetric neutrino sector must automatically contain the hidden Z| symmetry which 
dictates the solar angle 6 S as in 115. 21)) . 

As pointed out in Ref. [T], a particular choice of k = gives the conventional tri-bimaximal 
ansatz [56] tan# s = -7=- (9 S ~ 35.3°), but other choices of the group parameter k allow deviations 
from the conventional tri-bimaximal mixing, e.g., we can make a very simple choice of k = 3 , 
leading to tan# s = I (Q s ~ 33.7°), which agrees to the neutrino data equally well (cf. Table-1 
in Sec. 2) or even better (cf. Table-2 in "Note Added in Proof). The 7L 2 itself, as the minimal 
hidden symmetry for 6 S , is not restrictive enough to fix its group parameter k . But, extending 
the Z^ r ® Z| symmetry into a larger simple group can fix a particular k value and thus the solar 
angle 9 S . As we demonstrated in Sec. 6.3 of Ref. [1], a simple example is to enlarge Z^ T <S> Z| to the 
permutation group S4 [57], under which we can infer k = , corresponding to the tri-bimaximal 
mixing 6 S = arctan . 

Then, we examine how such a Z| symmetry could possibly survive after including general \i— r 
breaking terms in V' = V s + 8V . Expanding the small fj, — r breaking parameters up to NLO, we 
can derive the symmetry-transformation matrix G s corresponding to d v of (|5.12D . 

= G 2 = F'4 2 V't = V s dWv* + {V a d$8V* + 8V'd^ s ) 
= G° s + 5G S , 



G s 



(5.22) 



where 8G* 



Re[5G s 



Re[5G s ] + ilm[8G s ] with 
/ 

s 2 8 a +2c 2 s cos 5 D 8 X 



V2 

s 2a 8 a +2c 2 s cos S D 8 X 
V2 



s 2a 5 a +2c 2 a cos 5n&x 
V2 

2s 2 s 8 a -s 2s cos 8 D 8 X 




s 2s 8 a +2c 2 cos 8 D 8 X 
V2 



2s 2 8 a + s 2s cos 8 D 8 X J 



(5.23a) 
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/ n s 2s (8a-i—5a2)-2c 2 s sin5 D 6 x s 2s (Sa\ — 5az)+1c 2 s sinSpSx \ 

U V2 V2 

lm[5G s ] _ s ? j Sai -5a 2 y_2c* S in5 D 6 x q s *(Sa 2 - 8a 3 ) + s 2s sin 8 D 5 X 

y_ s? J Sai -S a ^ sSi nS D 5 x _ s 2 (6a2 _ 5a3) _ S2sSin5D5x o J 

(5.23b) 

where s 2s = sin 2# s . Because the symmetry transformation G s E Z| , we have the condition 
= X3 . Then, expanding this up to the NLO, we have verified the consistency condition, 

{G° s , 8G S } = 0. (5.24) 

Requiring that the Z| symmetry persists under fi— r breaking, i.e., the form of G s remains unaffected 
by the \i — t violation, we have the condition, 

G s = G° s , or, 5G S = 0. (5.25) 

Thus, with ()5.23ap -( |5.23b|) . we can derive the following solutions, 

= — cot 9 S cos 5d , (5.26) 

Ox 

from the real part condition Re[<5G s ] = , and 

25a l = 5a 2 + 5a 3 , (5.27a) 
5a 2 — 5a 3 = — 2 cot S sinfe S x = 2tan5 D 5 a , (5.27b) 

from the imaginary part condition Im[£G s ] = , where in the last step of (|5.27bj) we have made 
use of ()5.26p for simplification. Note that the correlation (|5.26p precisely agrees to what derived 
from our soft breaking model in Eq. (4.12a) of Ref. [1]; but now it is re-derived by requiring that the 
7j 2 symmetry persists in the presence of general low energy fi—T breaking. In addition, the above 
Eq. (|5.27b|) also coincides with Eq. (4.12b) of Ref. [I]. As we will demonstrate in the next subsection, 
the Z| symmetry is independent of the soft [i — r breaking in the seesaw model of Ref. [1] . We note 
that in the current construction of common fi— r and CP breaking with seesaw mechanism (Sec. 2. 2), 
such a Z| symmetry is not fully respected, hence the correlation (|5.26p no longer holds and we 
have predicted a modified correlation |^775] ), which can be tested against \5.2b}) by the on-going and 
upcoming neutrino oscillation experiments. 

To summarize, as we have demonstrated above from general low energy reconstruction formula- 
tion, the transformations G^ T = G 1 and G s = G 2 in the jjl— t symmetric limit correspond to the 
discrete groups 7j 2 t <g> 7L 2 , which are equivalent to and originate from the generic symmetry Z2 <8> Z2 
in the neutrino mass-eigenbasis because they are connected by the similarity transformations via 
(|5.13p . The fj, — T symmetry "Z 2 T has been known before, and the hidden symmetry Z| (as the 
minimal group dictating the solar angle 9 S ) was revealed by Ref. |T] in the context of neutrino seesaw. 
In this work, we further find that requiring the symmetry Z| to persist in the presence of most 
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general \x— t breaking terms will predict a new correlation (|5.26p between the small fi— r breaking 
parameters (5 a , S x ) . As we will prove below, the Z| symmetry is respected by a class of soft fi—T 
breaking seesaw models in Ref. p], but is partially violated in the present [i— r breaking seesaw model 
(Sec. 2. 2). 

5.2.2. Z| Symmetry and Neutrino Mass-Matrix with General fx — T Breaking 

In this subsection, we directly analyze the generally reconstructed light neutrino mass-matrix M v 
under the hidden symmetry Z| and the determination of solar angle 9 S . The mass-matrix (|3.ip can 
be uniquely decomposed into the fi— r symmetric and anti-symmetric parts, 

M v = Af W + 6M^ a) , (5.28) 

with 

(A B s B s \ /0 5B a -5B a \ 

M« = C s D , 5AfW = <5C a , (5.29a) 

\ C s j \ -5C a J 

B s = \ (Bi + B 2 ) , C s = i (d + C 2 ) , (5.29b) 

5Sa = \{Bi- B 2 ) , ,5C a = l - (d - C 2 ) , (5.29c) 

where we generally allow m l m 2 m^ ^ . Then, from ()5.9|) . the invariance equation of M v under G s 
corresponds to 

G\{M^ + 8M^)G S = AfW + <5AfW , (5.30) 

which uniquely gives, 

GtM«G, = Af«, (5.31a) 
GlSM^G s = SMj?) . (5.31b) 

Note that two possibilities may exist: (i). The Z| symmetry is a full symmetry of the light neutrino 
mass-matrix M„ if 6o£/i A5.31a\) and Ii5.31b\) hold. (ii). The Z| symmetry is a partial symmetry of 
M„ if the /x— r anti-symmetric part breaks (j5.31b[) . 

is) 

We can prove that the Z| is always a symmetry of the \i— r symmetric part M,; and generally 
holds (|5.31ap . Substituting (|5. 16|) into (|3.4j) and noting that the decomposition (|5.28p is unique, we 
can reconstruct the //— r symmetric and anti-symmetric parts of M„, respectively, 

AfM = V?D„Vj, (5.32a) 
5M^ a) = V*D V 5V' ] + SV 1 "^^ + 5^'* 

= F s * + SV 1 *!)^ + 0(5f ) , (5.32b) 
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where Sj denotes all possible NLO parameters under consideration (such as 5 X , 5 a and y' , etc). This 

(s) 

shows that the \i — t symmetric part M{, is diagonalized by V s . Hence, the corresponding Z| 
transformation matrix is just , as given by (|5.19p . The G® must be the symmetry of and 
thus always holds the invariance equation (|5.31ap . This proves that the solar mixing angle 9 S (as 
contained in the rotation matrix V s and symmetry transformation matrix G®) is generally dictated 
by the Z| symmetry, independent of any specific model. 

On the other hand, the validity of (|5.31bp is highly nontrivial because the requirement of G s = G® 
[cf. 15.25]) ] does not generally hold under \i — t breaking, and it has to be checked case by case. As 
we will prove in Sec. 5. 3, the fi— r anti-symmetric part will break Z| in the current /jl— t blind 
seesaw (Sec. 2), while it preserves Z| in the soft fji—r breaking seesaw of Ref. pQ. 

Using the expression of G s [Eqs. (|5.19ap and (|5.25j> ] . we can derive the solution from (|5.31ap for 
the fi — T symmetric part, 

2y/2 B s 

ta2 "- = A - ( c. + D) - (5 - 33) 
and another solution from (|5.31bp for the fi— r anti-symmetric part, 

tan0 s = -V2-^, (5.34) 



which further leads to, 



2y/2SB a 6C a 

t&n26s =~ 8ci-2sm ■ (5 - 35) 



Hence, if the Z| would be a full symmetry of M u (including its \i—r breaking part), the two solutions 
(I5.33P and (|5.35p for the solar angle 9 S must be identical, leading to a nontrivial consistency condition, 

. , 2V2 B s . 2V2 SB a 5C a 

tan2 ^ = A-(C S + D) ^ ~ 5CI-28BI ' {5M) 

An explicit counter example to this condition will be given in Sec. 5. 3. 2. 

In the following, we apply the most general reconstruction formalism (Sec. 3. 1) to compute the 
fjL—T symmetric and anti-symmetric parts of light neutrino mass-matrix M u = + M^ 1 ■ With 
these, we will explicitly verify Eq. (|5.33p by using the elements of fi — r symmetric , and we 

further derive physical consequences of the consistency condition (|5.36p by using the elements of fi— r 
anti-symmetric Mj) a \ 

♦Reconstruction Analysis for General Normal Mass-Ordering Scheme 



Eq. (|3.5p reconstructs all the elements of M v in terms of three mass-eigenvalues, three mixing angles 
and relevant CP-phases. The normal mass-ordering (NMO) has the spectrum m 1 < m 2 <C m 3 , so 
we can define the small ratios, 



m 9 m q — m 30 



yi = ^, y 2 = ^, y 3 ^^ ^- (5.37) 



m 3 
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Thus we have the independent NLO parameters for the NMO analysis, (y\, y%, z, 5 a , S x , Sa^ Scfrj) . 
Expanding them perturbatively, we derive the LO form of the fx— r symmetric mass-matrix M v , 

I \ 

= m 30 \ ~\ Y (5.38) 

with a 1Q = a 20 = a 30 = a Q , a 30 + 4> 3Q = nir , and the NLO elements in 5M U , 

SA = e" l2a o (e- i2 ^o c 2 yi + e- a ^s 2 s y 2 ) m 30 , (5.39a) 

5B ° = ^e- i2a ^e- i2 ^ yi - e -^2 y^ s[n2 e s m 30 , (5.39b) 

5C S + 5D = e~ i2a ° (e- i2 ^s 2 syi + e^fa^ m30 , (5.39c) 

SB a = -^e- iS n5 x m 30 , (5.39d) 

5C a = i [25 a — i (5a 2 — Sa 3 )] m 30 . (5.39e) 

From (|5.38p . we have ^4^ = = c[ ,-\-D^ = . Thus, using the //— r symmetric NLO elements 
(|5.39a|) - (|5.39c|) . we can compute the ratio, 

2a/2^ s {e- i2 ^ yi - e- i2 ^y 2 )^26 s 

A S -(C S +D) "" {e- i2 ^ yi -e- a ^y 2 ){ci-sl) " " 1 ' 

which explicitly verifies our Eq. (|5.33p [as generally derived from the invariance equation (|5.31ap 

under Z|] for the current NMO scheme. This is an explicit proof up to NLO that for a general NMO 

scheme the [i—t symmetric mass-matrix = Mj, + SMjf^ does hold the Z| symmetry. 

Then, using the \x — t anti-symmetric elements (|5.39dj) - (|5.39cp . we derive the ratio, 

-V2°-^ = °JR = tan^, (5.41) 

oC a 5 a — |(5a 2 — 5a 3 ) 

where in the last step we have used Eq. (|5.34p under the assumption that Z| symmetry also holds 
for the ijl — t anti-symmetric mass-matrix , i.e., the validity of the invariance equation (|5.31bp . 
Analyzing the real and imaginary parts of ()5.4ip . we deduce two relations, 

5 a = — S x cot 9 S cos So , (5.42a) 
Sa 2 — 5a 3 = 2tan5r)5 a - (5.42b) 

These are in perfect agreement with (|5.26p and (|5.27bp . which are generally derived under a single 
assumption that the Z| symmetry persists in the presence of ijl — t breaking. But, as will be shown 
in Sec. 5. 3. 2, this assumption does not generally hold, and the current [i— r blind seesaw (Sec. 2. 2) 
provides a nontrivial counter example. 

♦Reconstruction Analysis for General Inverted Mass-Ordering Scheme 

For the inverted mass-ordering (IMO), the light neutrinos have the spectrum m 2 > m 1 3> m 3 , so 
we can define the small ratios, 

z 2 = z 3 = — 5.43 
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where we have z x = z and z 2 — z + \y' in connection to the NLO parameters (y', z) introduced 
in Eqs. (|3.7|) - (|3.9|) of Sec. 3. 2. Then we have the independent NLO parameters for the IMO analy- 
sis, (z 1 , z 2 , z 3 , 6 a , S x , 8a^ <50j). Expanding them perturbatively, we derive the LO form of the 
symmetric mass-matrix M v , 



M<F> = mJ \ \ , (5.44) 




with a w = a 20 = a 30 = ot , </> 10 = (j) 20 = — a , and the NLO elements of M u , 

5A = m Q \(? s z l + s 2 z 2 — i2{c 2 5cj) l + s 2 5(j) 2 + (faq)] , (5.45a) 

5Bs = 27f m o sin 2e s[z!-z 2 -i2 (<50! -Sfa)], (5.45b) 

5C S + 5D = m [s 2 s z 1 + c 2 s z 2 - i (2s 2 5(j) l + 2c 2 5(p 2 + 5a 2 + 5a 3 )] , (5.45c) 

*£« = -^m e iS °5 x , (5.45d) 

<5C a = -m [S a + § (<5a 2 - Sa 3 )] . (5.45e) 

From (gSD , we have = and - (d 0) +D^) = 0. So using the ^-r symmetric NLO 
elements (|5.45ap -( j5.45cp . we can compute the ratio, 

2V2 B s sin2fl s [z 1 -z 2 -i2{5<t> 1 -5<t> 2 )} 

A s ~{C S + D) " cos 26 s \z x -z 2 -i2 {5<t> x - 5<j) 2 )} - i {2Sa x - 5a 2 - 5a 3 ) ' 1 ' j 

from which we deduce the consistent solution, 

a,% b :+d) - ^ 2e - < 5 - 47a) 

2<5oq = 5a 2 + 5a 3 , (5.47b) 

which explicitly verifies our Eq. (|5.33p [as generally derived from the invariance equation (|5.31aD 
under Z 2 ] for the current IMO scheme. Also the above solution ()5.47bp exactly coincide with the 
general Eq. (|5.27p . The above is an explicit proof up to NLO that for a general IMO scheme the fi— r 
symmetric mass-matrix Mv = + SJVfjf' does hold the Z| symmetry. 

Then, with the /i— r anti-symmetric elements ()5.45dp - (|5.45ep . we further evaluate the ratio, 

-V2°-^L = = tanfl s , (5.48) 

oC a 5 a + \ (5a 2 - 5a 3 ) 

where in the last step we have applied (|5.34p under the assumption that the /i— r anti-symmetric mass- 
matrix M{, also respects the Z| symmetry, i.e., the invariance equation (|5.31bj) holds. Inspecting 
the real and imaginary parts of (|5.48p . we deduce the following, 

5 a = — 5 X cot 9 S cos 5d , (5.49a) 
5a 2 — 5a 3 = 2tan(5cx5 a , (5.49b) 
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which coincide with Eqs. (|5.42ap -( |5.42bj) as we derived for the NMO scheme. We see that both 
(|5.49ap - (|5.49bp and (I5.42a|) - (l5.42bl) precisely agree to (|5.26|) and (|5.27bp which are generally derived 
under a single assumption that Z| is a symmetry of the full mass-matrix M v = M„ + 8M„ 
including its li—t breaking part 6M„ . But, as we will prove in Sec. 5. 3. 2, the above assumption is 
not generally true and for the fx— r blind seesaw with IMO (Sec. 2. 2) the Z| symmetry is violated 
by 5M ( „ a) . 

So far we have explicitly proven the relations (|5.26|) and (|5.27p for general NMO and IMO schemes 
via the general model-independent reconstruction formalism (Sec. 3.1), where the only assumption 
is that the Z^ symmetry fully persists in the presence of ll — r breaking. In the next subsection, 
we will map the Z 2 T <8> Z| symmetry into the neutrino seesaw Lagrangian, and demonstrate that 
the hidden Z| symmetry is a full symmetry of our soft li — t breaking model in Ref. [1] where the 
physical prediction (|5.26p holds; while for the current li—t blind seesaw model the Z| is only a partial 
symmetry (respected by the li—t symmetric part M„ ), and is violated by the li—t anti-symmetric 
part 5Ml a \ leading to our prediction of the modified new correlation (|4.13|) in Sec. 4. 1, in contrast 
to or Qg^gp . 

5.3. Mapping Z2 ® Z2 Hidden Symmetry into Neutrino Seesaw 

Consider the general seesaw Lagrangian in the form of (|2,2p with two or three right-handed neutrinos. 
After spontaneous electroweak symmetry breaking, consider the invariance of (|2.2p under the residual 
symmetry transformations, 

v L Gj -,u L , N -> GfN , (5.50) 

where Gj is 3-dimensional unitary matrix, and G^ is 2 x 2 or 3 x 3 matrix (depending on two or 
three right-handed neutrinos invoked in the neutrino seesaw). Accordingly, we have the following 
invariance equations for the Dirac and Major ana neutrino mass-matrices, 

GjrrijjGj 1 = m D , (5.51a) 
Gf T M R Gf = M R , (5.51b) 

from which we deduce the invariance equation for the seesaw mass- matrix of light neutrinos, 

G T jM v Gj = M v , (5.52) 

where M v = mj^M^m^ . Let us diagonalize the Majorana mass-matrices M v and M R as follows, 

VZM„U V = D u , UlM R U R = D R , (5.53) 

in which D v = diag(m 1 , m 2 , m 3 ) and D R = diag(Mi, • • •, M n ) with n = 2 for the minimal seesaw 
or n = 3 for three-neutrino-seesaw. Thus, from (|5.5ip - (|5.53p . we can express Gj and G^ as, 

Gj = U v d^Ut, Gf = U R d^U R , (5.54) 
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where kernel representation {d„ } is given in ()5.12p . and corresponds to the product group TL^ 

via the equivalent flavor representation {Gj} for the light neutrino sector. For d R in (|5.54p . we 

give its nontrivial forms, 

d!ft = diag(— 1, 1) , (for minimal seesaw), (5.55a) 

d^ft = diag(l, 1, — 1) , dft* = diag(— 1, 1, 1) , (for 3-neutrino-seesaw), (5.55b) 

where d R forms a Z,' 2 symmetry for right-handed neutrinos in the minimal seesaw, and {d$, d^} 

form a product group ® Z^ for right-handed neutrinos in the three-neutrino-seesaw. The trivial 

(i) 

case with d R equal to unity matrix is not listed here which corresponds to the singlet representation 
G R = X . Since the low energy oscillation data do not directly enforce a ll^ T symmetry for heavy 
right-handed neutrinos, we find two possibilities when mapping the Z^ r to the seesaw sector: (i). 
the right-handed neutrinos have correspondence with the light neutrinos in each fermion family and 
transform simultaneously with the light neutrinos under the Z^ T to ensure the invariance equation 
(|5,51ap ; this means Z^ 41 " = Z^ 1- . (ii). the right-handed neutrinos are singlet of the usual Z^ 1- 
symmetry (called ll fi—r blind"), so the extra symmetry Z^ in the M sector is fully independent of 
the Z^ T for light neutrinos; this means that under Z^ T the invariance equation (|5.51ap has G\ € Z^ r 
for light neutrinos and G R = X for right-handed neutrinos. As generally shown in Sec. 5. 2, the Z| 
symmetry dictates the solar angle 9 S for light neutrinos. The extra group T! 2 in the right-handed 
neutrino sector also has two possibilities: one is U 2 S = 7L\ , and another is for the right-handed 
neutrinos being singlet of the Z| symmetry with Gf = X . 

5.3.1. Neutrino Seesaw with Common Soft \l — t and CP Breaking 

In Ref. pQ, we studied the common soft fi— r and CP breaking in the minimal neutrino seesaw, where 
the right-handed neutrinos M = (Nu, N T ) T obeying the same Z^ 7 "(= Z 2 A "~) at the LO, and small 
soft fi — t breaking is uniquely constructed in M R at the NLO. In the fi — r symmetric limit, we 
inferred that the diagonalization matrix U R is a 2 x 2 orthogonal rotation with its rotation angle 
Or = 9^ = § [I], as expected. Thus, inputting ()5.55ap for d R , we deduce from (|5.54p . 

G* T = \ ; ~ I , (5.56) 

which is just the Z^ r transformation matrix for right-handed neutrinos. With the two right-handed 
neutrinos N = (Nu, N T ) T shown above, there is no rotation angle 9^ 2 and also no corresponding 

7J 2 symmetry. So the right-handed neutrinos can only belong to the singlet representation Gf = X 2 

(2) 

under Z| symmetry, with d R = X2 . In our soft fi — T breaking model [T], the Dirac mass-matrix, 

(a a\ 

m D = b c , (5.57) 





50 



exhibits the exact Z^ 1- symmetry, so it should obey the hidden Z| as well, 

G T S m D G* = m D , (5.58) 

where G s = G® is given by (|5,19p and G^ = X 2 ■ This further leads to the invariance equation for 
the seesaw mass-matrix of light neutrinos, 

G T s M y G s = M v , (5.59) 

where M v = m D M R 1 rn^ ) , and the invariance equation for Mr is trivial here since G^ = X 2 ■ [Given 
the form of G s = G® as constructed in (|5.19p . we can also explicitly verify the equations (|5.58p and 
(I5.59p .] Hence, the group parameter k of Z| and the corresponding solar angle 9 S via Eq. (|5.2ip are 
fully fixed by the elements of the //— r symmetric m D , and is independent of the soft \i—t breaking 
in M R (which is the Z| singlet). This is a general proof based on group theory, without relying on 
making any expansion of the /i— r breaking terms in M R . As can be explicitly solved from Eq. (|5,58p 
above, we have [T], 

tan# s = \k\ = tt^L (5.60) 
\b + c| 

As another nontrivial check, we inspect the consistency condition (|5.36p . With the form of M v 
in Ref. [I], we explicitly verify that ()5.36p indeed holds, 

2^2 B s 2y/2 5B a 5C a 2 v / 2 a (6 + c) 

tan2 ^ " A-(C S +D) ~ " 5CI - 25BI ~ 2a?-{b + cf ' (5 ' blj 

where both the \i—t symmetric mass-matrix M^f 1 and the anti-symmetric part SM^ determine the 
same solar angle 6 S . The last equality in (|5.6ip can be derived also from the solution (|5.60p above, 
they are all consistent. Hence, the Z| is a full symmetry of the seesaw sector and the light neutrino 
mass-matrix M v in this soft // — r breaking model. 

We note that this Z| symmetry has a nice geometric interpretation. The two vectors, u 1 = 
(a, 6, c) T and u 2 = (a, c, b) T , in the Dirac mass-matrix m D = (u 1 , u 2 ) , determine a plane S, 
obeying the plane-equation, 

x -^{y + z) = 0, (5.62) 

where the parameter k is given in (|5.2ip . As shown in Ref. [I], the 3-dimensional representation G s 
is just the reflection transformation respect to the plane S. For the case of three-neutrino-seesaw, 
the n — t symmetric Dirac mass is extended to a 3 x 3 matrix, 

fa' a a \ 

m' D = b' b c = (n , n l5 u 2 ) . (5.63) 
\V c b J 

Thus, to hold m' D invariant under the Z| symmetry, we just need to require its first column u = 
(a', b', b') T to lie in the S plane, i.e., 

-4- = ^ = k, (5.64) 
y/2V b + c K J 
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where k = tan# s as in ()5.2ip . This means that the Dirac mass matrix (|5.63p only contains one 
more independent parameter than that of the minimal seesaw; furthermore, m' D is rank-2 and thus 
detMj, = (det m' D ) 2 (det Mr) -1 = always holds, as in the minimal seesaw. 

5.3.2. fj, — t Blind Seesaw with Common fJ>—T and CP Breaking 

As constructed in Sec. 2, the fi — r blind seesaw defines the right-handed neutrinos N as singlet 

of Z 2 r symmetry. This means that we must have the Z 2 r transformation matrix G R T = Z 2 and 

(2) 

d R T = d R = I2 ■ Consider the general Dirac and Majorana mass-matrices in the minimal seesaw, 

fa a'\ 

Mr = I 

Mi 2 M 22 
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I) 



~ /M11 Mi 2 \ 

M R = L, J • ( 5 - 65 ) 



\6 2 c 2 J 

The Majorana mass-matrix Mr can be diagonalized by the unitary rotation Ur , 

UIMrUr = Mr = diag(Mi,M 2 ), (5.66) 

Then we can derive the seesaw mass-matrix for light neutrinos, 

M v ~ m D M^ l fh T D = mDM^mp , (5.67) 

where m D = rh D UR takes the form as in (|2.17p . For the jjl— t blind seesaw with N being Z^ 1 " singlet, 
we can always start with the mass-eigenbasis of N with M R = diag(Mi, M2) , which means that 
the rotation U R becomes automatically diagonal and real, Ur = X 2 ■ Then, the extra symmetry 
Z 2 of Mr must be independent of the of light neutrinos, i.e., Z 2 ^ Z 2 r . So the natural 
choice is Z 2 = Z 2 . The Z 2 can have a nontrivial d s R = d R = diag(— 1, 1) as in ()5.55ap . Thus, the 
corresponding symmetry transformation for M R is 

Gf = U R d s R U R = 4 1} = diag(-l, 1). (5.68) 

There is also a singlet representation of Z' 2 , corresponding to d R = I 2 ■ 

Then, let us inspect the possible Z 2 symmetry for the Dirac mass-matrix by including the fi — r 
breaking effects [cf. (|2.17|) and (|2.18p in Sec. 2.2]. This means to hold the invariance equations in 

G T s m D Gf = m D , G^MrG* = Mr, (5.69) 

which will become, in the mass-eigenbasis of right-handed neutrinos, 

G T s m D d s R = m D , d R T MRd R = Mr. (5.70) 

Since Mr and dR are both diagonal, the invariance equation for Mr always holds. So we can rewrite 
the above invariance equation for m D as, 

G T s m D d s R = m D , (5.71) 
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where m D = m D (m Mji) 2 . [The \i— r symmetric form of m D was given in Eq. (j2.8[) .] Using the 



notation m D , we can reexpress the seesaw mass-matrix, M v - 
deduce the invariance equations under G s and d s R , respectively, 



m (m n m 



D""D 



) . So we can further 



mr>mF D 



d iiu d u, r 



(5.72) 



Next, we inspect the two equations in (|5,72p to check the validity of the Z| symmetry after 
embedding the fx — r breaking into m D [such as those constructed in (|2.18|) for instance]. From 
(|5.72p , we will explicitly prove that the G s is a symmetry only for the fi — r symmetric part of 
M v oc (m D m^) ; while d s R is violated by the fi — r breaking terms in m^m D . Hence, the Z| 
symmetry is only a partial symmetry of the light neutrinos, valid for the \x—r symmetric part M u 

We can write down the mass-matrix m D with the most general fJ,—T breaking, 



m D 



/a 
b 



( 



+ 





8b 1 +8b 2 
2 

8b 1 +Sb 2 
2 



\ 



<5c 1 +<5c 2 
2 

<5c 1 +<5c 2 
2 



+ 





8b 1 —8b 2 
2 

y I Sb 1 -8b 2 



+ 



\ 

2 

(5"C]_— <5c 2 



— (°) _L X— W _|_ gfrffi) - — ( S ) _L X— (°) 



m" D ' + 5m D 



m" D ' + 5m D 



(5.73) 



where b 



8b l , b 2 



b-5b. 



2 ! 



5c 1 , and c 2 



5 c 



For the symmetric mass-matrix product, m D m D = M^/rriQ = M v , we compute, up to the NLO, 



/I 



m D m T D 





SO § 






(0 


a' a'\ 
2 2 


-(Sh+Sb 2 ) 


6 6 


- (Sci+dcz) 




C C 




V bj 






V 




/0 f -§> 






/o i - 


2 1 




b 


-(Sc t - 


5c 2 ) 


c 













\ 







i8b x -8b 2 



(5.74) 



where the M y denotes the sum of the first three matrices and <5Mj, equals the sum of the last 
two matrices. For deriving the LO matrix in (|5.74p we have used the relations (|2. 14j) for the 

IMO scheme. There exist two basic realizations for the common breaking of fi—r and CP symmetries 



in m D or m n : one is for 5b 



! = 5b 2 = and (<5c l5 5c 2 ) = c((' , Qe lu) ) , which corresponds to m D 
in (|2.18|) ; and another is for 5ci = 5c 2 = and (8b x , 5b 2 ) = (,e luJ ) , which corresponds to m D 
in ()2.23p . As we pointed out earlier, the invariance of the product (|5.74p under G s G 7L 2 [cf. (|5.72p ] 
would be justified so long as our general consistency condition (|5.36p could hold. So, with (|5.74p we 
can explicitly compute tan2# s from the two expressions in (|5.36p including the //— r symmetric and 
anti-symmetric mass-matrix elements, respectively. We thus arrive at 

2y/2B s a' 



tan 20^ 



A-(C S +D) 
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y/2t 



(5.75a) 



ta^> = -^jj = ^ = (5.75b) 



for 56i = = 0, and 



tan2tf s - §C 2_ 2SB 2 ~ a 2_ 2b 2 - tan4y S 



(5.76b) 



for <5ci = 5c2 = 0. The above explicitly demonstrates the inequality &i 7^ , and thus proves the 
violation of the consistency condition (|5,36p . This is because the \x — t anti-symmetric mass-matrix 
5M v a ^ = rh SM„ in (|5.74p breaks the TL S 2 symmetry. Hence, Z| is not a full symmetry of the 
mass-matrix M v . Nevertheless, we find that the fi — r symmetric part = fh Mj in (|5.74p 

does respect the Z| symmetry, and its invariance equation ()5.31a|) leads to the correct solution (|5.33p 

and thus (|5.75ap for the solar angle 8 S . Substituting (|5.'2ip into (|5.75a|) or (|5.76a|) . we derive the 

9 2 a ci 

equation, k -\ k— 1 = 0, with r = ^_ corresponding to t|5.75ajl . or r = corresponding 

to (|5.76ap . So we can fix the Z| group parameter k in terms of the ratio of seesaw mass-parameters 



k = -1± + (5.77) 
Finally, we compute the other symmetric product mJ,m D , up to the NLO, 

m T D m D = P J- {5b x + 5b 2 ) ( J " + 5c 2 ) r Y (5.78) 
\0 1 J \c 0J \b 2cJ 

The last two matrices of (|5.78p arise from the /jl—t breaking, which make m^m^, non-diagonal at 
the NLO, and thus explicitly violate the second invariance equation of (|5.72|) . This violation of Z| 
does not directly lead to observable effect at low energies since the seesaw mass-matrix M v for light 
neutrinos is given by the first product m mp in Eq. (|5.74p . Also, we could choose to assign the 
right-handed neutrinos to be singlet under the Z| from the light neutrinos, i.e., d R = X 2 , then the 
invariance equation for m D m r £ becomes trivial. But the first invariance equation in (|5.72p under 
G s 6 Z| is still broken by the {J.—T anti-symmetric mass-matrix 8M V °^ = m 5M^ in (|5.74|) for 
light neutrinos, as shown by Eq. (|5.75j) or (|5.76p above. 

From the analyses above, we conclude that the hidden symmetry Z| is a partial symmetry of 

(s) 

the present model, respected by the \i — r symmetric part M v of the light neutrino mass-matrix, 
and thus determines the solar angle 8 S as in Eqs. (|5.75ap and (|5.77p . This also agrees to the result 
p.ISp [Sec. 2. 1] or (j4.5|) [Sec. 4. 1] which we derived earlier. As a final remark, we stress that the 
violation of the hidden Z| symmetry by the (jl—t anti-symmetric mass-matrix SM^ = rh Q 5M„ in 



(|5.74p has an important physical impact: it predicts a modified new correlation and can be 



experimentally distinguished from Eq. 115.26]) as predicted before by our soft /jl—t breaking of neutrino 
seesaw [1] . 
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6. Conclusion 

In this work, we have studied the common origin of fJL—T breaking and CP violations in the neutrino 
seesaw with right-handed Majorana neutrinos being fi— r blind. The oscillation data strongly support 
jj, — t symmetry as a good approximate symmetry in the light neutrino sector, leading to the zeroth 
order pattern, (023, #13) = (45°, 0°) . Hence the fx — r breakings, together with the associated CP 
violations, are generically small. For the \l— r blind seesaw, we have convincingly formulated their 
common origin into Dirac mass matrix m D (Sec. 2. 2), leading to the unique inverted mass-ordering 
(IMO) of light neutrinos and distinct neutrino phenomenology. This is parallel to our previous 
work pQ where the common origin of fi— t and CP breaking arises from the Majorana mass matrix of 
the singlet right-handed neutrinos and uniquely leads to the normal mass-ordering (NMO) of light 
neutrinos. 

In Sec. 3, we gave the model-independent reconstruction of low energy fj, — r and CP breakings 
with inverted neutrino mass-spectrum. With this we derived various predictions of the [i — r blind 
neutrino seesaw in Sec. 4. In particular, we deduced a modified new correlation (|4.13p between the 
two small /j, — t breaking observables 623 — 45° and #13 — 0° , as depicted in Fig. 2 and is very 
different from that in Ref. pQ. Eq. f|4. 13j) is shown to also hold for the general three- neutrino seesaw 
in Sec. 4. 3. This correlation can be experimentally tested against Eq. (j4. 15aj) as deduced from our 
soft fji—r breaking seesaw mechanism [T]. As shown in Fig. 2 and Fig. 6, our predicted range of #13 can 
saturate its present experimental upper bound. Imposing the current upper limit on #13, we derived 
a restrictive range of the deviation, -4° ^ 6*23 - 45° ^ 4° at 90%C.L., in Eq. In Sec. 4. 2, 

we have further generated the observed matter-antimatter asymmetry (the baryon asymmetry) from 
thermal leptogenesis in the \i — t blind seesaw. Under the successful leptogenesis, we derived the 
constrained correlation between #23 — 45° and #13 — 0° , as presented in Fig. 6. This figure predicts 
a lower bound on the key mixing angle, #13 > 1° , which will be explored soon by the on-going 
reactor neutrino experiments at Daya Bay [TU], Double-Chooz [T2] and RENO [T3]. Fig. 7(a) further 
constrains the Jarlskog invariant J into the negative range, —0.037 < J < —0.0035 , while Fig. 7(b) 
predicts the range of neutrinoless /3/3-decay observable, 45.5 meV < M ee < 50.7 meV, which can be 
probed by the on-going neutrinoless /3/3-decay experiments [3] . A lower bound on the leptogenesis 
scale Mi is inferred from Fig. 4, Mi > 3.5 x 10 13 GeV, and is given in Eq. (|4.40p . The correlations 
of the leptogenesis scale Mi with the reactor angle #13 and the Jarlskog invariant J are analyzed in 
Fig. 9(a)-(b). 

Finally, we have studied the determination of solar mixing angle 612 and its connection to a hidden 
flavor symmetry Z| and its possible breaking in Sec. 5. The general model-independent Z2 <8> Z2 
symmetry structure of light neutrino sector was analyzed in Sec. 5.2.1. We first reconstructed the 
3-dimensional representation G® for Z| group in the /jl — t symmetric limit as in Eq. (|5.19|) . We 
proved that hidden symmetry Z| holds for any \i — t symmetric mass-matrix M v of light neutrinos 
and determines the solar angle 612 via its group parameter, k = tan 612 , as in Eq. (|5,2ip . Then, 
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requiring that Z| persists in the presence of general \i — t breaking, i.e., G s = as in ()5.25p . we 
deduce a unique correlation equation (|5.26p which strikingly coincides with Eq. (|4.15ap . as predicted 
by our soft fi— r breaking seesaw [I]. In Sec. 5.2.2, we further analyzed the validity of Z| symmetry 
from general model-independent reconstructions of light neutrino mass-matrix M u . We derived the 
general consistency condition (|5,36p for the validity of Z| symmetry in the presence of all possible 
fx — t breakings. Under this condition, we derived the nontrivial correlation (|5.42ap or (|5.49ap between 
the two fjt — T breaking observables #23 — 45° and #13 — 0° , which agrees to Eq. (|5.26p as derived 
earlier from pure group theory approach. We stress that the agreement between (|5.26p [or (|5.42ap ] 
and the prediction (|4,15ap from our soft fi— r breaking seesaw is not a coincidence. As we explained 
in Sec. 5.3.1, the true reason lies in the fact that the soft \i — r breaking is uniquely embedded in 
the right-handed Majorana mass- matrix Mr which is a singlet of the Z| group and thus does not 
violate Z| . On the other hand, for the fM — T blind seesaw, the fi— r breaking is solely confined in 
the Dirac mass-matrix m D which would have nontrivial transformation (|5.70p or (|5.72p if TL\ could 
actually hold. As we have verified in Sec. 5.3.2, the invariance equation (|5.72p hold only for the 
fi—T symmetric part of the light neutrino mass-matrix M u , and is partially violated by its fi—r anti- 
symmetric part [cf. Eq. (|5.74p ]. In consequence, we found: (i) the solar mixing angle 812 is dictated by 

(s) 

the group parameter k of the hidden symmetry Z| acting on the fi— r symmetric mass-matrix M u 
[cf. Eqs. (|5.75ap and (|5.77p ]; (ii) the consistency condition (|5.36j) no longer holds, and we predicted 
a modified new correlation ^.lSfy , which can be experimentally distinguished from Eq. (|4.15ap as 
predicted by our soft \x — r breaking seesaw [1] . In contrast to our previous prediction (|4.15ap , Fig. 6 
points to an important feature of the new correlation (|4.13p by showing a more rapid increase of 
#13 as a function of #23 — 45° [cf. also (|4.2Ua|) ]; this allows #13 to saturate the current experimental 
upper limit, and confines the deviation 623 — 45° into a more restrictive range, —4° ^ 623 — 45° ^ 4° 
at 90%C.L., as in Eq. (|4.22p . These distinctive predictions of the present \i—r blind seesaw can be 
systematically tested against those of our previous soft \i—t breaking seesaw pQ, by the on-going and 
upcoming neutrino experiments. 

Note Added in Proof : 

After the submission of this paper to arXiv:1104.2654 on April 14, 2011, two long-baseline accelerator 
experiments newly announced evidences for #13 via the v — > v e appearance channel, one by the T2K 
Collaboration [58J on June 14, 2011 and another by the Minos Collaboration |59J on June 24, 2011. 
Minos reported 62 e-like events above an estimated background of 49 events, and favors a nonzero 
#13 at 1.5(j level. The resultant confidence interval yields, sin 2 2Q\3 < 0.12(0.19) at 90%C.L. 
for NMO(IMO) with 5 D = 0; and the best-fit value is sin 2 20 i3 = 0.04(0.08) for NMO(IMO). 
On the other hand, the T2K experiment observed 6 e-like events with an estimated background of 
1.5 events, indicating a nonzero #13 at 2.5u level. This gives the 90%C.L. limits, 0.03(0.04) < 
sin 2 26>i3 < 0.28 (0.34) for NMO (IMO) with 5 D = 0; and the best-fit value is sin 2 20 13 = 0.11 (0.14) 
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for NMO (IMO). These new data indicate a relatively large #13 mixing angle, 



T2K: 



Minos: 



5.0° < 013(9.7°) < 16.0°, (for NMO), 

5.8° < 013(11.0°) < 17.8°, (for IMO); 

0° < 013(5.8°) < 10.1°, (for NMO), 

0° 013(8.2°) < 12.9°, (for IMO); 



(6.1a) 
(6.1b) 

(6.1c) 
(6.1d) 



at 90%C.L., where the central values are shown in the parentheses. We would like to point out 
that the new data from T2K and Minos further support our theory predictions which give the unique 
inverted mass-ordering (IMO) and favors a naturally larger 0i3 even for a rather small deviation of 
23 - 45° , as shown in Eq. ([4~2TIaD and our Fig.0 (Sec. 4.1) or Fig.© (Sec. 4.2). 

Shortly afterwards, a new global analysis of oscillation data has been performed [60] to include 
the latest T2K and Minos data. With this we can update our Table- 1 accordingly, and translate the 
improvements [60j into the new Table 2. 



Parameters 


Best Fit 


90% C 


.L. 


99%C.L. 


lcr Limits 


3cr Limits 


Am^lO^eV 2 ) 


7.58 


7.15- 


7.94 


7.07-8.09 


7.32-7.80 


6.99-8.18 


Am? 3 (lCr 3 eV 2 ) 


2.35 


2.20- 


2.55 


2.10-2.63 


2.26-2.47 


2.06-2.67 


#12 


33.6° 


32.0° - 


35.4° 


31.0° - 36.4° 


32.6° - 34.7° 


30.6° - 36.8° 


#23 


40.4° 


37.5° - 


47.9° 


36.3° - 51.3° 


38.6° - 45.0° 


35.7° - 53.1° 


#13 


8.3° 


5.09° - 


10.4° 


3.5° - 11.6° 


6.5° - 9.6° 


1.8° - 12.1° 



Table 2: The updated global analysis [BU] by including the latest data from Minos |59j and T2K 
long-baseline accelerator experiments. (Using the new reactor fluxes will slightly shift the mixing 
angles 0i2 and 0i3 a bit as shown in [60J.) 



With Table-2, we have systematically updated our numerical analyses in Sec. 4. We find that 
the predictions of Fig. [21 Fig. [6] and Fig. U] exhibit more constrained parameter space in an interesting 
way, while the other figures remain largely the same as before. For comparison, we use the updates 
in Table-2 and replot Figs.[21[6l[I] as new Fig.[10l Fig. [11] and Fig. [121 respectively. In Figs. [10] and [TT1 
we see that the updated 90% C.L. constraint on 13 (yellow area) just picks up the central region of 
our predicted theory parameter-space. Comparing the two plots in Figs. llOtfTTl we see that imposing 
successful leptogensis in Fig. II II makes the parameter space more centered along the two wings, and 
the region around 023 ~ 45° is clearly disfavored. Since the new global fit of Table-2 gives the 
90% C.L. limits, -7.5° < 23 - 45° < 2.9°, with a central value 23 - 45° = -4.6°, it is clear 
that the left-wing of the theory parameter-space is more favored over the right-wing. Furthermore, 
imposing the 023 and 0i3 limits from Table-2 on our parameter-space in Fig. 11, we deduce the 
allowed range at 90% C.L., —4.8° < 023 — 45° < 2.9° , which is shifted towards the negative side by 
about 1° as compared to Eq. (|4.22p . 
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Figure 10: Update of Fig. 2 (Sec. 4.1) by using the improved global fit in Table-2, with 2000 samples. 
The shaded region (yellow) shows the updated constraint on 8 13 at 90%C.L. 
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Figure 11: Update of Fig. 6 (Sec. 4.2) by using the improved global fit in Table-2, with 2000 samples. 
The shaded region (yellow) shows the updated constraint on 9^ at 90%C.L. 
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Figure 12: Update of Fig. 4 (Sec. 4.2) by using the improved global fit in Table-2, with 1200 samples. 

Then, Fig. [12] shows that the predicted parameter region in the M\ — 5 D plane is much more 
centered along the two edges in Fig. 4, and a high leptogenesis scale M\ > 10 15 GeV is strongly 
excluded except for the tiny regions of the CP-angle 5 D very close to 180° and 360°. 




Figure 13: Correlation of #13 and #23 — 45° as predicted by our Eq. (I5.26P [Sec. 5.2] without seesaw 
and under the assumption of exact Z| symmetry. Plot- (a) shows the correlation by using the global 
fit in Table-1, while plot-(b) depicts the correlation under the improved global fit in Table-2, with 
2000 samples in each plot. The shaded regions (yellow) give the allowed 90%C.L. ranges by the 
corresponding global fit. 

Finally, as a comparison, we further analyze the prediction from our Eq. (|5.26|) [Sec. 5.2] which 
we derived for light neutrinos alone (without invoking seesaw) and under the assumption of an exact 
Z| symmetry. As we already proved in Sec. 5.3, this Z| only holds in a class of models including our 
soft t and CP breaking model in Ref. [1], but can be violated in other class of models including 
the current fj, — T blind seesaw model. Hence, the Z| symmetry cannot generally hold in a model- 
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independent way. With Eq. (|5.26p . we plot the correlation between the \i — t breaking parameters 
#13 and #23 — 45° in Fig,ll3fa)-(b). In plot-(a) we show the correlation by using the global fit in 
Table-1, while in plot-(b) we depict the correlation under the improved global fit in Table-2. Each 
plots contains 2000 samples. The shaded regions (yellow) display the allowed 90% C.L. parameter 
space by the corresponding global fit. Note that Eq. (|5.26p holds for both normal mass ordering and 
inverted mass ordering of light neutrinos. Fig. 13 shows that our predicted parameter space can easily 
saturate the current upper limit on #13, and thus accommodates a relatively large #13 as indicated 
by the new data from T2K [58] and Minos [59] . The prediction of Fig. [13] differs from the above 
Fig. 10 significantly, because the coefficient in Eq. (|4.18ap [corresponding to Fig. 10] has a nontrivial 
suppression factor relative to that of Eq. (|4.18bj) or Eq. (|5.26p [corresponding to Fig. [13]. Furthermore, 
we note that the above Fig. [13] should also be compared to our previous Fig. 2 in Ref. [1] because the 
correlation (|5.26|) applies to both of them. But there are large differences between these two figures, 
the major reason is that we input the parameter #23 — 45° in Fig. 13 according to the oscillation data 
(Table-1 or Table-2) and without invoking seesaw, while the #23 — 45° in the Fig. 2 of Ref. [1] was 
derived as a function of fundamental \i— r and CP breaking parameters in the seesaw Lagrangian 
which were scanned within their theoretically allowed ranges. This also leads to a stronger upper 
limit of 6»i3 < 6° in Ref. []]. 

Note Added-2 : 

After the publication of this paper in Phys. Rev. D84 (2011) 033009, Daya Bay and RENO col- 
laborations announced new measurements of nonzero l3 on March 8, 2012 [61] and April 8, 2012 
[62], respectively. Daya Bay experiment made a 5.2cr discovery of nonzero # 13 |61j . sin 2 2# 13 = 
0.092 ± 0.016(stat) ± 0.005(syst); and RENO found a nonzero 6 l3 at 4.9cr level [62], sin 2 2# 13 = 
0.113 ± 0.013(stat) ± 0.019(syst) . These give the following 3a ranges of nonzero # 13 , 

Daya Bay: 5.7° < 6 l3 (8.8°) < 11.1° , (6.2a) 
RENO: 5.9° < 9 13 (9.8°) < 12.6° , (6.2b) 

where the numbers in the parentheses 6 l3 = 8.8° and 6 13 = 9.8° correspond to the central values. 

Then, we can re-plot the Fig.[TUl Fig.[TT] and Fig.fTSTb) as the new Fig.HH Fig.[TS] and Fig.fTUl 
respectively. In these new plots, we have scanned the experimental inputs within 3a ranges. For 
the successful leptogenesis, we find that the lower bound on the leptogenesis scale Mi becomes, 
Mi > 2 x 10 13 , at 3a level. The successful leptogenesis in Fig. [15] further requires, 9 13 > 1° . 

To compare with our predictions, we have displayed the 3a range of 9 13 from the new Daya Bay 
measurement [ST] in the green shaded region. Furthermore, we show the 3a lower and upper limits 
of #13 from the new RENO data [62] by the horizontal red-lines. The horizontal black dashed-lines 
in each plot denote the 3a limits from the global fit |60j . From Figs. ll4lfT5l we see that the new limits 
from Daya Bay |61j and RENO [62] experiments nicely pick up the central regions of our predicted 
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parameter space of 6 13 . 




023-45° 

Figure 14: Update of Fig. 10, with 2000 samples. The experimental inputs are scanned within 3a 
ranges. The 3a ranges of the new Daya Bay data |61| are shown as the green shaded region; and the 
3a limits of the new RENO data [62 j are depicted by the horizontal red- lines. The horizontal black 
dashed-lines denote the 3a limits of the global fit [60J . 




Figure 15: Update of Fig. 11, with 2000 samples. The experimental inputs are scanned within 3a 
ranges. The 3a ranges of the new Daya Bay data |61| are shown as the green shaded region; and the 
3a limits of the new RENO data [62J are depicted by the horizontal red- lines. The horizontal black 
dashed-lines denote the 3a limits of the global fit [60J . 
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Finally, Fig. 1 161 shows the correlation of 8 13 and 623 — 45° as predicted by our Eq. ()5.26p [Sec. 5.2] 
under the assumption of an exact Z| symmetry and without invoking seesaw. We have scanned the 
3a ranges of 6 12 in (|5.26p . We see that the new data of Daya Bay [61] and RENO [62] pick up the 
upper parts of our predicted parameter space of 9 13 . 




023-45° 

Figure 16: Correlation of #13 and 623 — 45° as predicted by our Eq. ()5.26p without seesaw and under 
the assumption of exact Z| symmetry, with 2000 samples. The experimental inputs are scanned 
within 3a ranges. The 3a ranges of the new Daya Bay data [61] are shown as the green shaded 
region; and the 3a limits of the new RENO data [62j arc depicted by the horizontal red-lines. The 
horizontal black dashed-lines denote the 3a limits of the global fit [60] . 
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